arXiv:hep-ph/9502241v3 9 Feb 1995 


IC/95/15 

|hep-ph/9502^ 


An Introduction to the 
Heavy Quark Effective Theory 

F. Hussain and G. Thompson 
International Centre for Theoretical Physics, Trieste, Italy 


Abstract 

These lecture notes begin with a brief survey of the physics of heavy 
quark systems. This discusion motivates the introduction of the Heavy 
Quark Effective theory (HQET) which captures a great deal of the intu¬ 
ition developed. A derivation of the HQET from QCD is presented as 
well as an analysis of its special properties. The effective theory can be 
seen to amount to a one-dimensional field theory in the quark sector. The 
heavy quark flavour- and spin- symmetry of the effective lagrangian is an 
offspring of this. Other topics covered include the question of covariance 
of the theory, the construction of interpolating fields for the heavy hadron 
states, the application of LSZ reduction theorems to determine the (re¬ 
duced) number of form factors in flavour changing transitions, a complete 
verification of Luke’s theorem, plus the matching conditions between QCD 
and the HQET beyond tree level. 

These are an expanded version of lectures presented during the Trieste 
1994 Summer school on High Energy physics. 
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1 Introduction 

QCD as the theory of strong interactions has been with us for over twenty years. It 
has been remarkably successful in describing high energy physics. The discovery 
of asymptotic freedom has allowed for many perturbative calculations of physical 
quantities within QCD when combined with the parton model, such as Drell-Yan 
processes, deep inelastic scattering and various scaling relations just to name a 
few. These successes notwithstanding QCD has proved to be intractable in the 
infrared regime by direct analytic means. For example we do not have a handle on 
conhnement. One, at this point, normally needs to employ lattice gauge theory 
or to pass to QCD inspired effective theories such as non-linear sigma models 
or potential models of various kinds. Put rather blankly-QCD in this regime is 
hard. 

However, during the last few years there has been a resurging interest in 
heavy quark physics within the context of QCD. The reason for this is that 
one has been able to extract general principles from particular models of heavy 
flavour transitions. In the work of |l| one can see the importance played by the 
velocity of the heavy quarks (as opposed to their momemtum); indeed hidden in 
the analysis is a peaking about a velocity super selection rule It is to the 
credit of Isgur and Wise [Q that they were able to extract the physical, model 
independent implication, namely, that for inhnitely heavy quarks, velocity is the 
only important parameter. 

Further, 

(i) there is a heavy quark flavour symmetry, 

(ii) there is a two-fold spin degeneracy (because the spin coupling is oc l/mg, 
which tends to zero), and 

(iii) at the zero recoil point, or equivalently at maximum momentum transfer, 
the elastic transition is absolutely normalized (following Voloshin and Shifman 

i)- 

Our understanding of the dynamics has not improved. However, the new sym¬ 
metries in heavy quark theory give rise to numerous predictions for free! 

The heavy quark effective theory HQET, as it has come to be known, cap¬ 
tures the physics of the heavy quark systems which brings to light these new 
symmetries. HQET, as we will see, is an expansion of the QCD action in inverse 
powers of the heavy quark mass. Heavy quarks Q, for present purposes will be 
(c, b, t) with masses mq ^ (1.5, 5,175) Gev. The scale is set by Aqcd — 330 Mev 
and the expansion parameter is typically Aqco/'mQ- See Table 1. While the top 
quark seems to be a theorists delight, the expansion parameter being about 10“^, 
it unfortunately decays far too fast to be useful. 

Working within the context of lattice theory Eichten and Hill had suggested 
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spin 0 spin 1 

niD = 1-86 Gev niD* = 2.01 Gev 

= cu or = cd D*^ = cu or D*^ = cd 
niB = 5.28 Gev uib* = 5.32 Gev 

B~ = bu or = bd B*~ = bu or B*^ = bd 


Table 1: Gharm and bottom mesons. 


a non-covariant effective action for the infinitely heavy quark limit of QGD, based 
on previous work of the related, but somewhat different, non-relativistic limit of 
QGD by Feinberg ^ and by Gaswell and Lepage |^. This action was covari- 
antised by Georgi ^ and subsequently a great deal of effort has been expended 
on determining the implications of this theory, both on the lattice and in the con- 

H, 


tmuum. 


One of the important developments has been the observation 


by Luke pT[ that in the case of heavy meson transitions the Voloshin-Shifman 


normalization at zero recoil remains valid even when l/mq corrections to the 
lowest order effective theory are taken into account. This no-correction theorem 
has been generalised and proven in various ways in the literature |^, ||22||. 


A derivation of the effective theory which also gives a systematic expansion 
to any order in l/mg was given in |^, and Previously, corrections to 
the inhnitely heavy quark effective theory had been obtained by matching the 
results to QGD. As in the heavy quark effective theory is obtained from QGD 
directly by a series of Foldy-Wouthuysen type of held redehnitions, the matching, 
at tree level is immediate. The same is true for the derivation given in [Q. That 
one needs to ‘improve’ the relationship between QGD and the HQET at one 
and higher loop level comes about because regularization does not respect the 
derivations. 

These notes are based on an unpublished manuscript written with J. G. Ko- 
erner and S. Balk and on a talk delivered at the 1992 DESY workshop given 
by the second author. As the reader will see, we have also benehted from the 

m 


review articles of M. Neubert 


and B. Grinstein 


The notes have turned 


out rather long on account of the fact that we attempt to give a rather com¬ 
plete account of the properties of the heavy quark effective theory at 0{1) and 
0{l/mQ). It has been our intention to give a critical account of the assumptions 
that go into the formulation of the HQET. We have tried to be clear on what is 
derived and what is intuited, pinpointing the assumptions that are being fed into 
the analysis. Our hrst aim then is to make precise what we can, with conhdence, 
say about heavy quark physics. The following picture emerges. 


1.1 Physics of a Heavy Hadron 

The classical picture of a heavy hadron that one has is in many ways similar to 
that of the Hydrogen atom. The typical momentum carried by the light degrees 
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of freedom (light quarks plus glue) inside a hadron Aqcd is of the order of the 
proton mass divided by three ~ 330Mev. This is a measurement of how far 
the light quarks are “off” shell or, rather more correctly, this is about what 
their ‘constituent’ masses are. Typically, then, the light quarks and their gluonic 
cloud are carrying momentnm Aqcd- In heavy hadrons the ‘mass’ of the heavy 
qnark uiq is much greater than the typical scale, uiq 3> Aqqd, and the heavy 
quark is carrying most of the momentnm of the heavy hadron. The interactions 
of the heavy qnark with the light degrees of freedom will also only change the 
momentum of the heavy quark by the order of Aqcd, so that the heavy quark is 
then almost on mass shell. Indeed one expects Mq ^ uiq + O(Aqc'd)- 

While the change in momentum of the heavy quark is of order Aqcd, its 
change in velocity, Aqcd/^q -C 1, is negligible as the mass of the heavy quark 
goes to inhnity. The picture one has then is of the heavy quark moving with 
constant velocity and this velocity is that of the heavy hadron. It is important to 
emphasise that it is not momentum that is being eqnated but, rather, velocity. 
In the rest frame of the heavy hadron the heavy quark is almost at rest; it is only 
slightly recoiling from the emission and absorption of soft glnons. This picture 
does not depend on the actual value of uiq but just that it satishes uiq 3> Aqcd- 
As the mass of the heavy quark is taken to be bigger and bigger the recoil is less 
and less nntil ultimately, in the limit tuq oo, the heavy qnark does not recoil 
at all from the emission and absorption of soft gluons. In this limit the heavy 
qnark acts, therefore, like a static colour source. It is clear that in this limit the 
binding is independent of the flavour and hence the difference between the mass 
of the heavy hadron and the heavy quark, A = Mq — tuq, is a universal, flavonr 
independent, constant. 

In many ways we have also just decribed the Hydrogen atom. Take the proton 
to be a fnndamental particle. The typical momentum imparted to the proton by 
the electron and the photon cloud is very much smaller than the mass of the 
proton {nip). The proton can be taken to be a static photon source which binds 
the electron to form the Hydrogen atom. The fact that the Schrddinger equation 
for the electron in a 1/r potential describes qnantitatively the Hydrogen atom so 
well is indicative of the snccess of this pictnre. One of the things that is missed 
by this ‘non-relativistic’ analysis is the hyperhne splitting of energy levels. But 
such corrections are rather small compared to the enery level, AE/E 1. One 
can derive the Schrddinger equation from the fully relativistic and interacting 
Dirac theory for the Hydrogen atom and systematically incorporate corrections 
snch as the Thomas term for the spin-orbit coupling. 

Likewise for the heavy hadron, it is immaterial, in a hrst approximation, what 
the spin state of the heavy qnark is and, in analogy to the above discussion, one 
can give a systematic derivation of corrections to this pictnre. The corrections 
will be of the form of a series in l/niQ with the spin coupling coming in at next to 
leading order. One of our aims in these notes is to give a derivation from QCD of 
this ‘effective’ theory whose lowest term is analogous to the static proton action 


4 



in QED. 

We do not need, however, to get into any highly abstract formal considerations 
in order to exhibit the features mentioned above within the context of QCD 
1^ . Let us simply take the rriQ —>■ oo limit of the lowest order propagator and 


connected three point function for the heavy quark near its ‘mass-shell’, that 
is with pq = ttiqv + k, where v is the velocity of the heavy hadron and the 
components of k are bounded by Aqcd- The hrst of these behaves as. 


^Q-rriQ 




( 1 ) 


This asymptotic form already exhibits some of the features of the physical ar¬ 
guments. Recall that QCD with massless quarks has a flavour symmetry. This 
symmetry is (softly) broken by the introduction of different masses for the differ¬ 
ent flavours. Clearly there is no dependence on the mass at lowest order in (|^) 
and so, just as for massless quarks, there is no dependence, at this order, on the 
flavour of the quark. 

The second important feature is that the lowest order term is diagonal in 
spin. The (1 -|- '/) projector picks out the particle state in the held while (1 — ff) 
projects out the anti-quark component of the held. The rest of the propagator is 
proportional to the unit matrix and so does not depend on the spin of either the 
quark or anti-quark components. 

The three point Greens function (we amputate the gluon leg) 


i 

^Q-rriQ 




i 

/q + 


with soft gluon momentum q (also bounded by Aqcd) goes like 


( 2 ) 


. + OiAgcD/mg) . (3) 

2 v.kv.{k + q) 

This expression also picks out the quark component and does not depend on spin 
to leading order. 

To summarise: the heavy quark symmetry is the statement that 


HEAVY HADRONi(spini) = HEAVY HADRON 2 (spin 2 ) 

1.2 Analogies and Differences 

We can get some more information as regards the physics of heavy quark theory 
by pushing the analogy with the Hydrogen atom further. At some point the 
analogy must breakdown and we will have to come to grips with that as well. 
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Recall that the Schrodinger equation makes no reference at all to the nucleon 
and so it is invariant under the symmetry described above. For the moment let us 
think of the different isotopes of Hydrogen (p+e), Deuterium (pn+e) and Tritium 
(pnn+e), as arising from different ‘flavours’ of the proton. As far as the electron 
is concerned the photonic held is the same regardless of the havour. Treating the 
nuclei as inhnitely massive the Schrodinger equation describes these systems well. 
The spectra of the three (ignoring reduced mass corrections) are identical. This 
coincides with the heavy quark symmetry which allows one to exchange heavy 
quarks without effecting the spectra. Of course the Schrodinger equation is blind 
to the spin of the electron relative to the proton. At this order of approximation 
this is the other part of the heavy quark symmetry, namely that the spin of the 
heavy quark is immaterial. 

What the Schrodinger equation misses is the hyperfine splitting of the energy 
levels. A more correct account follows from an analysis of the Dirac equation. 
Still treating the nuclei as very heavy the predictions of the Dirac equation differ 
from those of the Schrodinger equation at order l/rrip. The correction of interest 
is a Pauli-term with the typical coupling. The corrections can be 

systematically derived by a series of Foldy-Wouthuysen transformations. We will 
apply these to QCD. 

Now it is time to mention the major difference between QCD and QED. 
Field theoretically we know that ag (e^/dvr) and ag {g'^/ATi) run. For large dis¬ 
tances (Ue tends to zero while ag grows. The reason that one can perform with 
some reliability the calculation of the spectra of the Hydrogen atoms is that 
the coupling constant is small. Now the typical distance is the Bohr radius 
a = /i^/(27re)^me ~ .5 x = 5 x 10^ fermi so the binding occurs at “large” 

distances compared to the typical length scales of the nucleons which are of the 
order of fermi’s. However, in QCD although, again, binding is a long distance 
phenomena, the coupling constant is large and perturbation theory is therefore 
not reliable. 

Recall that, in the static proton approximation, the potential for an electron 
in the photonic held of the proton depends on distance as 1/r. One obtains this 
behaviour in the Born approximation by consideration of hgnre 1. 

The corresponding amplitude is 

up{p + k)'ygUp{p)Df^^{k)Ue{q - k)'y„Ue{q) (4) 

where p' — p = k, and p = mpv and v = (1, 0). Apply now the Gordon decompo¬ 
sition, 

up{p + k)Yup{p) =up{p + k) i d- up{p) (5) 

y Zuip IfTlp J 

on the proton side to arrive at 

up{p + k) ( ^ Up{p)Dg„{k)Ug{q - k)Yug{q). (6) 

y 2mp 2mp J 
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In a covariant gange the photon propagator is 

. ( 7 ) 

Substituting this into (H), and remembering that u{p + k)^u{p) = 0, one obtains 

up{p + k)'y^up{p)-^Ueiq “ k)'y^Ueiq) ■ (§) 

Now, as p^ = p'"^ = rup, we find that 2mpko + k^ — k'^ = 0. In the infinite 
mass limit, so as not to have an infinite momentum transfer, one finds that 
kf) ^ k‘^/2mp. Notice that in the heavy proton limit, k/mp 0, the Gordon 
decomposition goes to up{p + k)v^up{p), which is consistent with our tree level 
analysis of the vertex in that limit. As /cq —^ 0, (||) reduces to 

- k)j%e{q) + 0{l/mp). (9) 

—k^ 

This reproduces the usual result of the coulomb scattering of an electron in a given 
external electromagnetic field (generated by a heavy nucleus) A^{x) = {ie/A-nr, 0). 

It is only the 0-component of the gauge field which enters finally so that one is 
tempted, as we will be tempted in due course, to ask what happens in the gauge 
Ao = 0? Physics is gauge invariant, so that going to this gauge changes nothing. 
Indeed the photon propagator in the gauge v.A = 0 is 

D^y{k,V)-{p^y (10) 

and when sandwiched between the currents only the p^y part survives. The 
calculation is the same as in the covariant gauge and one reproduces (P). 

A problem arises when one passes to the heavy proton limit first. The naive 
limit allows us to replace with 

up{p-\-k)v^up{p) (11) 

but, in the gauge v.A = 0, the photon propagator satisfies v^D^y{k,v) = 0, so 
that at first sight, it seems that there will be no scattering at all! The resolution 
of this puzzle is that in non-covariant gauges spurious poles in the propagator 
mix terms of different orders in 1/mp. Let us first dispense with the term up{p + 
k)afj,yk'^Up{p)D^'"{k,v)ue{q — k)paUe{q) in (H). By current conservation and anti¬ 
symmetry of a^y the only part of the propagator that survives is that proportional 
to Pfj,y. This is the same as for the covariant gauges and can safely be ignored in 
the limit. In passing to the limit (|lTl) we have dropped a term 

up{p + k)^^up{p) ( 12 ) 
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which is clearly Xjrtip down. However, 




k'^ — v'^v.k 

{v.ky 


(13) 


We can disgard the k^ piece by the usual argument, to hud 

2mpv'^ 


k^D^^{k,v)^- 


v.k 




(14) 


While (p!^ is indeed 1/mp down it may not be neglected as the photon propagator 
is not just of order 1 but there are components of it that are of order mp. 

The upshot of this discussion is that one cannot simply apply the heavy quark 
limit blindly. We have seen that while one is, at the level that we are working, 
allowed to use the covariant gauges one cannot employ the axial gauge without 
mixing different orders in the heavy quark (proton) mass. 


1.3 Physics of Flavour Changing Transitions 

We will be interested, in these notes, in calculating the flavour changing electro- 
weak transitions of one heavy hadron into another. The transitions of prime 
interest will be where the heavy quark in the hrst interacts with the electro-weak 
particle to flavour change into the heavy quark of the second. The light quarks 
will essentially be spectators. The reason for this is that in the heavy mass limit 
the heavy quark, as we have seen, is essentially on-shell and acts as a colour 
source for the light degrees of freedom. In particular the spin of the heavy quark 
decouples from the dynamics. Thus the dynamics of the weak transition of a 
heavy hadron are essentially determined by the point like interaction of the weak 
current with the heavy quark. The picture which emerges of the transition is as 
follows. 

In the inhnite mass limit in the rest frame of the heavy hadron, the heavy 
quark Q is at rest and is surrounded by the light cloud with no spin interaction 
between them. In the transition the heavy quark Q emits a W meson and becomes 
another heavy quark Q' moving with some velocity (the velocity of the hnal heavy 
hadron). The light cloud has thus to adjust its velocity to keep up with the heavy 
quark Q' in forming the new hadron. It is this adjustment, or overlap, of the light 
degrees of freedom which gives rise to the form factors. We immediately see from 
this simple, atomic physics, picture the physical reason why the unique form 
factor in heavy meson decays, for example, is normalised to one at zero recoil (or 
q^ax)- this kinematical point, the daughter heavy quark is produced at rest 
in the initial rest frame. It is clear that nothing has changed as far as the light 
degrees of freedom are concerned because there is no flavour dependence in the 
static colour force due to the heavy quark and the light degrees of freedom do not 
feel the effect of the change in the heavy quark mass (both are inhnite) unless it 





moves. Thus there is a complete overlap of the light wavefunctions before and 
after the transition at this kinematical point and hence the form factor is one. 
Put bluntly there is no ‘dynamics’ in the transition at this point. 

2 Effective Theory to Lowest Order 

It is our aim in this section to derive an effective held theory from QCD which 
captures the essence of the heavy mass limit. In this limit the mass of the heavy 
hadron, the bound stateQ, is approximately the mass of its heavy constituents, 
but it is important to remember that these do not make up the total momentum 
of the heavy hadron; although taken to be nearly massless, the light degrees of 
freedom are allowed to carry off momentum. One is tempted to say that, if there 
is just one heavy quark inside the hadron, the 4-velocity of the bound state is 
the 4-velocity of its heavy quark. This is not quite true, but something can be 
said about the heavy quarks 4-momentum projected in the direction of flight of 
the heavy hadron, as will be made concrete below. 

2.1 Basics of the Effective Theory 

Georgi ^ describes the situation rather vividly. If the bound state, made up of a 
light quark system and one heavy quark, is moving with velocity = 1, > 

0), then the 4—momentum of the bound state is 

= MqV^^ = TTlQVf^ + F , (15) 

where Mq is the mass of the bound state, essentially the same as the mass 
of the heavy quark mq (as indicated by the second equality). Clearly = 
{Mq — mQ)v^, which in the weak binding limit is = T,qmgV^ where rriq is the 
light quark mass. If the bound state undergoes scattering into a new state of the 
same heavy quark, with momentum 

PS = mQv'‘‘+ k’‘‘, (16) 

(n'^ = 1), the momentum transfer is mQ{v — v')^ + {k — k')^. In the limit as 
mq —> oo, for hxed k^ — k'^, then necessarily v'^ = to keep the momentum 
transfer hnite. Conversely, to alter the velocity of the heavy hadron takes an 
inhnite amount of 4—momentum. So one is able to follow the velocity of the 
heavy hadron. Soft gluons will not alter the heavy bound state velocity. Only 
very hard gluons or electroweak transitions can do this. 

^Subsequently, unless otherwise stated, whenever we say the heavy particle (or heavy 
hadron) we mean the bound state containing the heavy quark. Mq is always the mass of 
the Q-hadron while mq is the Q quark mass. 
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To sharpen this picture, consider the heavy quark’s behaviour during this 
scattering event. The heavy quark is off shell, but as the physical mass is (in 
the limit tuq —> oo) identihable with the heavy quark mass, the heavy quark 
cannot be ‘far’ off shell. It carries most of the four momentum of the system. 
Decompose the hadron 4—momentum, before the scattering, as a sum of the 
heavy quark 4—momentum pq and that carried by the light degrees of freedom 
Pi 

Pq = mqv + k = pQ+Pq. (17) 

The heavy quark four momentum is then pq = mqv + [k — Pq). Repeating 
the argument above, now for hnite Pq as well, establishes that the heavy quark 
4-momentum after the scattering event is pq = rriqv + [k' — Pq) {p'q hnite). 

From this we learn that the velocity of the heavy quark is not hxed, but rather 
the projection of the heavy quark’s four momentum in the direction of the four 
momentum of the hadron is the heavy quark mass, as for hnite {k — pq), v.{k — pq) 
is very small compared to the mass. Thus 


v.pq ~ mq. 


( 18 ) 


The idea then, is to construct an ehective held theory that is able to keep track 
of heavy hadrons with given velocity and which also encodes the fact that the 
heavy quark momentum satishes (plSj) “on shell” . 

Georgi |^, following Eichten and Hill [Q, did precisely this in the limit that 
the heavy quark masses are inhnite. The action that he hnds in this case is 

- mq)Q , (19) 


for a given heavy quark with a given hadron velocity and with the covariant 
form of the projection ([T^) as the kernel^, where D is the covariant derivative 
— igA^. Here A = and T“ are the usual Hermitian SU(3) Lie 

algebra generators. The heavy quark effective theory at this order will be called 
HQETq. 

We have seen, in the introduction, that the heavy quark propagator can be 
expanded as 

^ .(l + ^)/2 


^Q-mq 


= 


v.k 


0(k/mQ). 


( 20 ) 


This is precisely what one gets for the zeroth order term on expanding, in powers 
of 0{k/mq), the tree-level heavy quark propagator 


i 

^v.pQ - niQ 


( 21 ) 


^ It is not quite this as we have not rotated out the mass, nor have we summed over all 
possible velocities, but the actions are essentially equivalent. 
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arising from the effective action (p!9D. In fact one can also write 

1(1 + :^) . 1 ( 1 - 1 ^) 


]/v.pQ - rriQ 


= i- 


— i- 


v.pQ - niQ v.pQ + tjiq 

. 1(1 + . 1(1 - 

= t- -^- t- 


v.k 


( 22 ) 


v.k + 2mQ 

and then one notes that the heavy quark mass pole is in the first term and as 
tuq —cx) the second term goes further and further away from the pole. All the 
higher order, in (l/mg), terms arise from the second factor in the propagator 
(P^. Subsequently in Green’s functions when considering terms of zeroth order 
we shall often take the heavy quark propagator to consist only of the first term 
in the above equation (p^. 


Similarly, the 3-point function (a heavy quark and a gluon) is given by 


Gf^\pQ.q) = 


^Q-rriQ 


i'igi’' 


^Q + 




(23) 


where A^^{q) is the gluon propagator. Now consider soft gluons, i.e. q is also of 
order Aqcd- Then we can expand as before to obtain 

.1 + 




\ r i^^/'l + 


V {k + g)'^''^(^) G{Aqcd/ mq) . (24) 


Exercise: Check the form of the two and three point functions in the heavy quark 
limit. Hint: 


1 + ^ 
2 


)7^( 


2 ^ 



1 + 
2 


(25) 


This corresponds precisely to the vertex, obtained from the effective 

Lagrangian ([T^), when sandwiched between the leading order term in the heavy 
quark propagator (|22|) . Thus in the limit mq oo the only part of the gluon 
held which contributes is that along the n-direction. These results can be easily 
extended to arbitrary tree level diagrams, provided that we have only one heavy 
quark and all other particles carry small momenta. 

The theory dehned by ([I^ has some rather remarkable properties. All of 
the simple characteristics and enlarged symmetries that the theory enjoys may 
be traced back to the fact that the operator {i]/v.D — mq) that appears in (^) 
depends only on a one-dimensional derivative and only on one linear combination 
of the Dirac matrices, namely For all intents and purposes this action dehnes 
a one dimensional held theory in the heavy quark sector (whence solvable in that 
sector). This should be kept in mind. 
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2.2 Flavour Symmetry 

The flavour symmetry at this order, valid in the static limit, may be seen as 
follows. The combined action for a b and c quark with hadron velocities Vb and 
Vc respectively is 

+ S°{vc) = j b{ifSb'Ub-D - mb)b + J c{ifl^Vc.D - rricjc . (26) 

The flavour mixing symmetry is between b and c quarks in hadrons moving with 
the same velocity, so that the relevant form of (|26|) is 

S^{v) + 5'°(u) = J b{i'j/v.D — mb)b + J cii'^v.D — mc)c , (27) 

where v = Vb = Vc- The symmetry transformation is 

6c = eexp {i]/v.x{mb — mc))b, 

6b = —eexp {i'^v.x^rric — mb))c, (28) 

where e is an infinitesimal parameter. If one includes the top quark (or other 
heavy flavour quarks) there is a corresponding set of transformations leaving the 
action invariant. 


Exercise: Derive the corresponding invariance for F heavy flavours. 


In order to be rid of the cumbersome exponentials in (^) one may, if one 
wishes, define new fields 

Q = exp{i'j/v.xmQ)Q , (29) 

in terms of which the action (^) has the same form but the mass is set to zero 

i 

S°q{v) = j Q{i,lv.D}Q. (30) 


This property that one may define fields in terms of which the mass is transformed 
away is characteristic of the one-dimensional first order operator. Eq.(^) is the 
infinitesimal form of a more general global symmetry, here a U(2) symmetry (in 
general U(F) for F heavy flavours). Thinking of the {b, c) as an U(2) doublet the 
action 


E 

Q=b,c'' 


Q{i]6v.D)Q 


(31) 


is invariant under 


(6,c) ^ {b,c)U 



(32) 


where U eU(2). 
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2.3 Particle Properties 

A second characteristic of the one dimensional nature of the held theory is the 
fact that one may create ‘covariant’ velocity projection operators to project out 
quark and anti-quark states. Let 

P±(t>) = i(l ± /), (33) 


with the properties 

Pi = P+, Pl = P-, P±P+ = 0 , P+ + P_ = 1. (34) 

Now decompose the heavy quark held as 

Q = Q+ + Q- , Q± = P±Q 5 (35) 

so that Q+ satishes P-Q+ = 0, which is the Dirac equation (/— mQ)Q+ = 0 for 
a quark with momentum tuqv^ and Q- satishes the Dirac equation for an anti¬ 
quark with the same momentum. As operators these annihilate precisely those 
states. Note that these projection operators act in much the same way as the 
parity projection operators do in the rest frame of a ‘normal’ particle, which when 
boosted become Dirac equations . In any case the action (|T^ decomposes into 
quark and anti-quark pieces 

Sq{v) = J Q+{iv.D - mQ)Q+ - Q_{iv.D + mQ)Q_ . (36) 

These projections are also reminiscent of chirality, and do not tell us about the 
dynamics of the wave equation itself (except of course that it has two orthogonal 
sets of solutions). The free wave equation for Q+, for example, is 

{iv.d - mQ)Q+= 0 , (37) 

with the solutions 

Q+aix) = P+„ , (38) 

where F+ is an arbitrary Dirac spinor depending only on = x^ — v^v.x which 
satishes x^.v = 0, so that F+j^^x-^) is unaffected by v.d. 

Physically, the situation is then somewhat underdetermined. Certainly the 
wave function satishes the Dirac equation (its ‘chirality’ is fulhlled) but what 
happens orthogonally to the line of hight of the hadron is not known[|. The 
reason for this is not too surprising. A one-dimensional held theory can only 

•^In previous works we have pointed out, rather euphemestically, that it is the heavy quark 
‘label’ that satisfies the Dirac equation with velocity v and not that the heavy quark itself has 
only that velocity. We hope our remarks here clarify the situation. 
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determine the dynamics in that direction. There is a similar solution for Q- with 
similar conclusions. 

An alternative way of saying this is the following. The fact that the wave 
equation has such arbitrariness in its solution is an indication that the action has 
a rather large invariance. We may as well exhibit this for (plQl). Let 


Q ^ M{x^)Q, Q ^ QM 


(39) 


with M diagonal in colour space, but a general matrix in spinor space. To 
guarantee invariance of the action one only needs that [y', M] = 0, which is 
partially solved by M = P+M_ + P_M_|_ where M_|_ and M_ are functions. These 
invariances of the action account for the arbitrariness in the wave functions (|3^ . 
A complete decomposition of the invariances is given in the appendix. 


2.4 Covariance and Spin (I) 


How covariant is the construction thus far? The actions (|T^ and (|36D depend 
explicitly on a hxed four-vector, so that arbitrary Lorentz boosts cannot be a 
symmetry of the action. Having chosen a preferred 4—velocity we wish to consider 
the subgroup of the Lorentz group which keeps this time-like direction hxed. In 
the heavy hadron’s rest frame = (1,0) the actions are clearly invariant under 
the 0(3) group of three space rotations. Let a general inhnitesimal Lorentz 
transformation be parameterised by 7]^^ -|- where is the metric tensor. 

Then the little group is dehned to be generated by all those which satisfy 




= 0 


(40) 


It is easy to see that, for satisfying this condition, the heavy quark action is 
invariant under the transformations 

SQx = , (41) 

with the consistency property that the particle projection operators P± commute 
with the spin transformation, [P±, = 0, so that the transformed helds 

remain unambiguously eigenstates of P±, i.e. 


P^5Q± = 0 . (42) 

The spin transformations (l4l|) require some explanation. There is no trans¬ 
formation (Lorentz rotation) of the co-ordinates x^. This comes about because 
the operator {i'^v.D — mg), as we noted before, acts only in the direction of v, 
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so any Lorentz transformation ^ yx'^ with ''Vy = 0 is by itself a 

symmetry of the action. For the helds this translates into 


5Q±{x) = yX''d^Q±{x), 

5Q±{x) = -X^ yx''d^Q±{x), (43) 

5Ap = -X^yx'^d^Ap (44) 


which is a symmetry of the action. The two sets of transformations (|4T|) and 
together form the conventional Lorentz transformations for the spinors. 


Exercise: Check the invariance of the heavy quark action (|^ under the above 
transformation rules. 


Note also that the transformations (|4l|) , the solutions given after (|39|), and 


= lb^^iv)MQQ , 5Q = (5/-^(r;)75Mo , (45) 

with Mq a function of x^ and ^^{v) = ^ — y'n.e ^ exhaust the symmetry of the 
action (|^) in four dimensions. Indeed, the transformation (^T]) may also be 
generalized as 


5Q = -Q|f^(n)75^Mi, (46) 

where 7 is dehned by X^'^'j^jy = and Mi is a function. 

2.5 Field Choices and Interactions 

As only one component of the gauge held appears in the action 5'Q(n), namely 
v.A, if we pick the gauge v.A = 0 then the heavy quarks associated with that 
hadron velocity decouple from the glue. Such non-covariant gauges may be quite 
difficult to work with, and in any case are of limited use for they provide no 
simplihcation for heavy particles with other velocities. In fact, as we saw in the 
introduction, this gauge can, in certain situations be dangerou^. We would like 
to see this decoupling, therefore, in a more general fashion. Just as it was possible 
to dehne helds for which the mass does not appear in the action (|^) , we can 

"‘For a general tensor of rank n, we set (r^) = {i’)Ayi....^y„ 

with T]^y{v) = - VpVy. 

® Though for properties of transitions at equal velocity, that is at maximum momentum 
transfer, such gauges have been very useful p^. 
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similarly define fields Q and Q for which the gauge interaction is absent. These 
are given by 


Q{x) 

Q{x) 


where the Wilson line 


W 


Q(x) 


Q(x)W 


-1 


W 


X 

V 


Pexp 


W 


dsA ■ V 


(47) 


is a path-ordered exponential, wherein the path is a straight line from cx) to a; 
along the n-direction. 

In terms of these fields ([l^ becomes 


Sg(v) = J QXilv.d - mo)Q . 


(48) 


So Sq represents free heavy quarks in any gauge. Let us explain how this comes 
about. To say we can pick a gauge v.A = 0 is to say that given an arbitrary 
gauge field A there exists a gauge transformation, which takes us to that gauge. 
Specifically, if 


ki' = r'^e + r'-5e, 

9 

(49) 

and v.A' = 0 then 


0 = -I- -v.d)^ , 

(50) 

9 

is an equation for solved by 


e(n.a;) = W 

(51) 


with the boundary condition that, at v.x = —oo, ^ is the identity (actually 
we could just as well take ([501) to be the defining equation of the path-ordered 
exponential). Substituting (^71) into (|T^) we find 


= J Q{i]^v.D-mQ)Q 
= J -mQ)^Q 

= J Q{i]/v.d - mQ)Q + i J - igv.A^)Q 

= J Q{i]/^v.d - mQ)Q . (52) 


We see that,in changing variables according to (^Tj) , the gauge field A is altered 
to the gauge field A! which satisfies the gauge condition v.J^ = 0, giving the 
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triviality of the action. So there is no gauge condition on A. It is just that 
the combination which makes up A' satishes (|50|) . The helds Q are gauge inert, 
meaning that they do not vary under gauge transformations, as can be seen from 
their dehnition. This ensures that the action (^) is gauge invariant as it should 
be. 

Let us note that this is rather an extreme situation. The effective theory is 
totally disinterested in the coupling of gluons to the heavy quarks. This shows 
us that the action (^) really models the situation described above, namely, that 
while keeping track of the heavy particle velocity, soft gluons are not important 
for questions concerning transitions. 

In terms of the effective helds (|47|) the action (|^ of HQETq takes the free 
form described above (|^. The symmetry ( ^31) is thus directly generalised to the 
free case by 

5Q{x) = yx'"d^,Q{x) , 5Q{x) = yx'"d^,Q{x). (53) 


For completeness, let us note that the Feynman propagator for the non¬ 
interacting helds Q is 


Sf{x - y]v) 


-i^^^-^0{v.x - v.y)5l{x - 

~ ^\ {v.y - v.x)6l{y - 


(54) 


where the argument of the delta function (5^(a;) is the component of x not in the 
n-direction. To pass to the Feynman propagator for the interacting held Q is 
particularly simple in view of (^71), 


SF{x,y;v, A) = Sf{x 


y;v}W 


W 



(55) 


Notice that these only propagate quarks forward in “time” (we mean v.x which 
is the usual time component in the rest frame) while only anti-quarks propagate 
backwards in “time”. One immediate consequence of this is that there are no 
heavy quark loops in the HQFTq, as could be inferred directly from the fact that 
they are certainly not there in the v.A = 0 gauge. 


Exercise: Show that in QED one can calculate any Green’s function involving 
Q and A^ non-perturbatively by going over to the free spinor variable Q. Hint: 
Calculate a Wilson loop in pure Maxwell theory. Why does this procedure not 
work in QCD ? 

2.6 Derivation of Lowest Order Action 

Questions that naturally arise are, how does one get to the effective action (^) 
directly from QCD and what are the corrections that are needed for hnite mg? 
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As recognised in |^, this action is familiar in a slightly different but related 
context. Namely it is a generalisation of the so-called ‘Bloch-Nordsieck’ model 
that one arrives at in studying the infrared behaviour of QED as the 
electron goes on shell. The similarities of the two contexts in which this action 
arises will not have escaped the readers attention. Indeed we know how to derive 
the non-relativistic action, and we may use this knowledge to get at the heavy 
quark effective theory. 

Recall that to obtain the relevant non-relativistic information about the Dirac 
equation one makes use of Foldy-Wouthuysen transformations. These are canon¬ 
ical unitary transformations which eliminate the interaction terms between the 
positive and negative frequency components of the fermion spinors in the Hamil¬ 
tonian. A peek at ( p!9tp6D , and the subsequent discussion, shows us that this is 
precisely the situation that we would like to get to. 

In held theory the analogous set of transformations that give a systematic 
derivation of the non-relativistic limit of the relativistic Hamiltonian in non- 
Abelian theories was presented by Feinberg P (see also |^, |^). The non- 
relativistic limit picks out a preferred frame (usually the rest frame, but this 
need not be so). Likewise, in the heavy quark limit, one would like to hx or pick 
out the “velocity” of the h- (or other heavy) quark that matches that of the b- (or 
heavy) hadron. Furthermore we would like to arrive at a systematic expansion (in 
the inverse mass) giving corrections to this picture. This is achieved by making 
use of the Foldy-Wouthuysen transformations in the context of Lagrangian held 
theory. 

Here we derive (|^) by ignoring l/mg corrections, which will be taken up 
again in later sections. Consider the action 




ip - mQ)V'Q, 


(56) 


where the subscript Q, as before, indicates a heavy quark. We make no distinction 
here between renormalized and unrenormalized quantities, as it is not relevant 
for our present purposes and would clutter up the formulae. Following we 
perform the change of variables 


= ^i{P-iv.D)/{2racPQ 

= - ^n.D)/(2mQ)Q + 0(l/my + ... , 

= Q-Qi(^5)/(2mQ)+0(l/mg)+ ... , (57) 


where D is dehned by 

j fDg^-j fDg. 


(58) 
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The Jacobian of this transformation within dimensional regularization has been 
shown to be unity in . We stay within the realm of dimensional regularization 
in the following. 

The action may now be expressed as: 


Sq = J Q[l + iiP - ]^v.D)/{2mQ)]{ip-mQ)[l + i{p - ]^v.D)/{2mQ)]Q 
= j Qiiiv.D - mQ)Q + Oil/rriQ) + ... 

= Sl){y)+ 0{l/mQ) + ... , (59) 


which is the relationship we have been looking for. A feature that has just been 
used and that will persist throughout is that the mass in (^) lowers the order of 
the terms in the transformation (^Tl), allowing for the replacement of the covariant 
derivative in (|^ ) with the one-dimensional covariant derivative. 

The reason for the particular exponential form of the Foldy-Wouthuysen trans¬ 
formation will be explained in detail when we tackle the derivation of the effective 
theory at 0{l/mQ). For now, we note that to really be able to denote the extra 
terms in (|5^) as Op/mq) we must be sure that there are no components of the 
derivatives that appear that he in the v direction. This is exactly what the form 
of the Foldy-Wouthuysen transformations (|57|) guarantees. 

In this derivation we have assumed that both, the components of the mo¬ 
mentum of the heavy quark, and the components of the gauge held that are 
orthogonal to the line of hight are small relative to the (inhnite) heavy quark 
mass. For this to be the case, the interactions must be such that there are no 
hard gluons involved. 


2.7 Covariance and Spin (II) 

While, as we have discussed above, the effective action (^) is invariant under 
little group transformations it is clearly not invariant under general Lorentz trans¬ 
formations. 


5xQ{x) = LxQ{x) 

= (-A'^ ^x^df, + ^X^''li,-fu)Q{x) , 

SxQ = LxQ 

JaA, = -X^.x'^d^Ap + X/A,. (60) 

This circumstance prompted Georgi [Q to integrate over all possible velocities, 
with (|T^) as the integrand, so as to restore manifest covariance. Somewhat later 
this integral over velocities was demoted to a sum although it was not clear what 
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the domain of the sum should be. The derivation, presented in the last section, 
of the effective theory directly from the QCD action, fixes once and for all the 
four-velocity to that of the physical hadron. 

The non-covariance of the action ([I^) under the Lorentz transformation (| 6 ^) 
is actually a desired feature. This may be seen from two points of view. Firstly, 
think of the Lorentz transformation in the active sense, that is, as a boost of the 
hadron from its old four velocity u to a new four velocity v' = u-fA.u. The effective 
action designed to describe the new situation is Sq(v'), which is precisely what 
one obtains under the transformation (|60|) , namely Sq(v') = S'q(v) — (5a5'q(u). 
There is no transformation on v in the variation. This was, after all, to be 
expected. 

What we really want are the usual covariance identities that arise from Lorentz 
invariance; for example (with 

(//e(«;)|lV„6](0)|i/,W)) = A„''(ff,(«,)|pJ„6](0)|i/,(«,)). (61) 

The velocity dependence on the right hand side rests in the action SqIvq). Per¬ 
form a Lorentz transformation on all the helds with parameter A^ j, on the right 
hand side. This is just a change of dummy variables in the path integral ap¬ 
proach. The states IHqIvq)} map to IHqIvq)}, while the action, as we saw, 
becomes Sq(vq). Thus (|6lD is established. 

Now to the second point of view. The effective action taken to all orders 
is equivalent to the QCD action, which is covariant. Hence covariance must be 
achieved in the full effective theory. But incorporating higher order (in I/uiq) 
terms in the action would hardly help the situation. The missing piece is that 
(|60|) is not the correct transformation rule for the effective helds! If P~^{v) = 
P — ytv.D transformed correctly as a bi-spinor then the effective helds would 
have the transformation rule (| 6 ^) . However, as u is a hxed vector, we hnd, using 
the dehning equation ( 0 . to hrst order in l/mg, that 


SxQ{x) 

S\Q 


LxQ{x) + -^[v.X.'yv.D '^v.X.D]Q{x ), 

- --^Q[v.X.'^v. D +^u.A. D] . 

2mQ 


(62) 


Here a.X.b = a^X^'^b^. It is straightforward to establish that (l 6 ^), alongwith 
the usual transformation of the gauge held, leaves the lowest order action (^) 
invariant upto order l/mq. Taking into account higher order terms in both the 
action and the transformation allows one to establish covariance at any given 
order in the inverse mass. 


Exercise: Derive (| 6 ^) . 
mations. 


For fun determine the l/rug corrections to the transfor- 
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3 Current Induced Transitions 


For definiteness consider the transition of a bonnd state containing one 6-qnark 
|$b) to another bound state containing one c-quark |<hc), induced by a flavour 
changing current J^ = V’c7At(l — 75 ) "06 = Then we can write 

the transition-matrix element as 




(63) 


The situation may be described as follows. As the b hadron comes in from the 
far past, we are able to track its four velocity Vf,, and after the transition we may 
follow the c-hadron’s four-velocity Vc into the future. In this way the relevant 


action to consider is (p6D 


S^{vb) + ^ / b{i^b''^b-D - mb)b + / c{i^^Vb.D - mb)c , 


(64) 


which is repeated here for convenience. 

Let us embroider on the picture of the heavy quark theory that has been 
developed thus far. A bound state made up of one heavy quark and some light 
quarks moving with velocity v, such that the 4-momentum of the system is 


= Mqv^ , 


(65) 


undergoes a current induced transition to a heavy bound state of a different 
flavour, with a new four velocity v' and four momentum 


= M, 


Q'V 


ifi 


( 66 ) 


Then the momentum transfer, = Mqv^ — Mqiv'^] in the limit as the heavy 
quark masses go to infinity it diverges. The square of the momentum transfer is 
= Qhax - ‘^MqMQiiv.v' - 1) and is bounded by = (Mg - Mqif, which is 
achieved at equal velocity v = v'. Of course, even though q has components that 

are large, may well be small e.g. = 0 at n = (1, 0) and v' = ( 2 m m^' ’ —) with 


V = 


2MqMq. 


)-,0,0) say. The velocities are well defined when the ratio Mq/Mqi 


is finite, which we take to be the case. 


3.1 LSZ Reduction and the Effective Theory 


To get a handle on the transition (|63D we first make use of the ideas of the 
reduction theorems. These tell us that as long as we can And an ‘interpolating 
field’ ((){x) corresponding to a state |*h(F)) with the property (for the F-wave 
pseudoscalar state) 

(O|0(a;)|$(F)) = e-*^-, (67) 
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( 68 ) 


then the transition may be expressed as^ 

where 

L(fi, Pc) = lim - M;)(P; - M^) . 

We shall concern ourselves with the exact form of the interpolating helds 0 q 
presently. For now we note that whatever their precise form, they will have the 
structure 


(t)Q{x) = Nq,QXp^{d)T'tpf{x)'tl)Qia{x), (69) 

(I>q{x) = Nq.QXa'^{d)T^Q{x)'ipqip{x), (70) 

with X = 7oxWo for heavy mesons, while for heavy baryons one has 

(f)Q{x) = iV$Qe*^^x“^^(a)TV’Qia(a:)^/’gj73(a;)^/>q/fc^(x), (71) 

(^gix) = N^QeijkXap^{d)Tiljll{x)ijf{x)ij’l^^{x), (72) 


with X = x^7o7o7o; where x°‘^"' are projection operators (Dirac matrices possibly 
including differentiations, total and partial, like d) which pick out the particle 
states of interest^. The normalization constant for brevity denoted by Nq, 
is hxed by (|^ . ifjg and ^|Jq are the heavy and light quark helds, respectively. 

To simplify the formulae we concentrate on mesonic transitions, i.e. from a 
Pmeson to a c-meson. The generalisation to baryons is immediate. One may 
rewrite (^) in this case as 

NbN^L J 

(0|W“(4)^g(2:)^/’ca(2:)^M(0)X/3^('9j;)^f(2/)V’gp(2/)|0) • (73) 

Now any total differentiations with respect to x and y, present in x{dy) and 
x{dx), can be shifted on to the exponentials by partial integration and, calling 
the subsequent projection operators xi^b) and xi^c), one can write the matrix 
element as 

NbNMPb:Pc) I d^xdSP^-^e-^P-y 
{^\^b{y)-x{vb)-'4^q{y)^q{x).x{vc).iJc{x)Jy,{^)\{^ , (74) 

®A11 the matrix elements are described in terms of path integrals. 

^Here, for once, we explicitly show the trace over the colour indices i,j,k. They will be 
implicit usually. 
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From now on (.) also indicates snm over Dirac indices, which should be clear 
from the context. Since we will nearly always consider interpolating fields inside 
reduction formulae with the relevant exponentials, we may as well always replace 
x(d) with y;(n) in the interpolating field from the beginning. From now on we 
shall do so unless otherwise stated. Thus far we have followed the full theory in 
conventional form. Now we apply the effective field theory. In order to do so, 
we need to assume, that in the complete theory, the integral in ([74|) will factor 
according to 


(P2 _ M2) (P2 _ M2) 


r/° + 


1 1 1 1 1 N 

+ 0{—^) + 


2m 


2mh 


m 


(75) 


Q 


where is the term that arises in the lowest order effective theory and tjI and rjl 
denote first order contributions. Let us naively apply the effective field theory. In 
(|74D replace everywhere the heavy quark fields ipQ with the corresponding fields 
in the effective theory Q{x) = W 


Q{x) to give 

N,NMPb,Pc) J d^xd%e^^-^e-^P-y 

{0\b{y).x{vb).M{y,x).x{vc).c{xW).T^.b{0)\0), 


L^-qJ 
4 


with 


X 0 0 y 

Vc Vc Vb Vb 


(76) 

(77) 


M{y,x)f, “ = Tri,,g{y)i,‘{x)W 

where the trace is over the colour indices. Here we have defined the matrix 
product 

W r^i^2^3a;4l ^ pp Ixil pp Ixal-lpp Ixal pp ra;4l"^ /'yg'j 

\_Vi V2 V2, V4, \ L ”^1 J L ”^2 J L ”^3 J L "^4 J ^ 

We have used the fact that the heavy quark propagators are diagonal in 
colour. Note that M{y,x) depends explicitly on Vb and Vc- As the effective fields 
Q decouple from the interactions, (0 factorises into 


NbNMPb^Pc) j d"a;dV^-"e-*^*-*'Trx(n),)(0|M(|/,a;)|0). 


x(n,)(0|c(a;)g(0)|0)F^(0|6(0)%)|0). 


(79) 


Set the momentum through {0\M{y,x)\0) to be q and k. Then the transition is 
easily expressed as 


NbN^L{Pb, Pc) j d^kd^qTrM{q, k) 
Xivc)-. - ^ .. -T,, 


cVc-{Pc -k) -mc^^ ib'^b-{Pb + q) -rrib 
where M[q, k) is defined through the Fourier transform 

{0\M{y,x)\0) = [ d^qd*ke-^'^-ye-^’^-^M{q,k), 


X{vb) , 


(80) 


(81) 
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and we have introduced the heavy quark propagator 

(0|Q(a:)Q(|/)|0) = / -. (82) 

J pv.p - niQ 

3.2 Taking the on-shell limit 

To ensure that the transition (|80D does not vanish there must be poles generated 
by the integrals that match the zeros that arise from the factors outside the 
integral. This is an assumption about bound states that needs to be made in 
the full theory before the effective theory is invoked. However, this need not 
be taken for granted within the context of the effective held theory, as we have 
mentioned around (|75|) ; rather, due to the simplihcations encountered we may 
and will address this issue in some detail. 

We are not free to take the limit in any fashion we choose. For example, let 
us set Pq = MqVq + ^q^vq), for some inhnitesimal eg, and understand the limit 
to be 


L{Pq) 





(Recall that e-*-^(n) = — v^v.e so that n.e-*- = 0.) The problem with taking the 

limit in this fashion is that Cg cannot appear anywhere in (|5(]| ) except in L{Pq)\ 
The mass pole is not regulated in this way. An immediate consequence is that 
we have no choice but to take the limit along the vq direction. We then take the 
limit in the direction Pq = Mg(l + eg)ng, so that 


lim (P^ — Mq) = lim 2egM^ 


(83) 


We can leave the cancellation of the zeros implicit in (|80|) , but nevertheless 
simplify its form. The explicit form of x(n), that will be given in the next 
section, satishes — l).x(n) = 0. The propagators as contracted against the 
wave functions then take the form 


-mq) y 


1 


VQ.{PQ-k) - mg 




(84) 


so that the transition is quite generally 

N,NMPb,Pc) f d^kd^q— 

J Vr- 


{Pc -k) -rUc Vb.{Pb + q) -rUb 


TrM{q, k)x{vc)Tf,x{vb) 


(85) 


There is, however, a way to take the limit so that the expected cancellation of 
the pole and zero may be explicitly achieved. This requires that the heavy quark 


24 







propagator {vQ.{PQ—k)—mQ) ^ exactly cancels the zero {Pq—Mq). This requires 
that the Green’s function (0|M(g, k) |0) is peaked, within the integral, at Vc-k = 
Ac and —Vb-q = Ab, in the heavy quark limit, where Aq = Mq — mq. Then, for 
example, the c-quark propagator becomes {vc-Pc — me — Ac)~^, exactly cancelling 
the zero [P^ — M^). These conditions are precisely the ‘on-shell’ assumptions that 
we made in determining the heavy quark theory in the hrst place. We assumed 
in (IB) that the projection of the quark momentum in the direction of the hadron 
4-velocity is the quark mass and in this limit the hadron mass. If this cancellation 
is accepted then (IS^) becomes 


iNtNMM, / d‘kd‘qTrM{q,k)x{v,)T.x(vi) 


( 86 ) 


It is not mandatory that the cancellation of the poles proceeds in this way. For 
our purposes it is enough to notice that the denominators of the heavy quark 
propagators depend only on Aq and on eqMq. However, we know (see Introduc¬ 
tion) that to lowest order Ab = Ac = A. So that after the limits are taken the 
only dependence will be on A, which is universal. 

A problem that still remains is, how do the light degrees of freedom know 
that there is such a momentum restriction (i.e. vq.k ~ A)? As we have noted 
the information that is lost in going to the effective fields has to do with how the 
light degrees of freedom respond to the heavy quark. The discusion thus far is 
then still not complete. We need to ask ourselves about the applicability of the 
effective theory in such a naive fashion. Certainly at the point of interaction with 
the current the only important scale is that set by the transition, and the heavy 
quark operators may be replaced with their effective counterparts. But that still 
leaves the heavy quark operators that go into making up the interpolating field. 

We parameterise this ignorance by the substitution 

Xaai...an ^ Xal3i...l3n^ai...ani^) ’ 

SO that the spinor functions B ‘correct’ the behaviour of the light quarks and soft 
gluons to remind them that they are in the presence of the heavy quark, leaving 
the heavy quark label a ‘free’. This choice corresponds to the physical picture 
presented in the Introduction that in the heavy mass limit the heavy quark spin is 
decoupled. These functions should be considered to contain the nonperturbative 
information (both in the coupling and in the mass), that is needed for a bound 
state to emerge. One should expect that, while these are certainly functions of 
the heavy mass, there may be no sense in which they can be expanded in a 1/mq 
series. For example, such a B might have a dependence of the form 

xjmq/fi^exp {-mq{vq.kf / , (88) 

which is a sum of derivatives of a delta function in a 1/mq expansion and which, 
consequently, may yield quite divergent integrals term by term. This is reminis¬ 
cent of the non-analytic behaviour one expects in terms of the strong coupling in 
the confining region in QCD. 
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If one accepts, for the moment, the need to introduce the S-functions into 
the dehnition of the wave function, then it is apparent that the discussion on the 
normalization of the wave function needs to take this into account. One would 
have to insert B in the appropriate places in (|134|) , (|135|) , (|136|) and (|138|) . 

With the introduction of the B factor, the form of the matrix element (^) 
does not change but now M (y, x) is dehned as 




(89) 


As M{x,y) is in any case unknown, this ‘improvement’ will not effect any of 
our subsequent discussion on the number of form factors. One diagonal contri¬ 
bution to the function B is the correction coming from the matching of HQET 
to QCD at one loop level for example (see Section 6). 


3.3 Bargmann-Wigner Wave Functions and Interpolating 
Fields 

In order to extract useful information from the overlap integral for the transition 
we need to determine an appropriate set of projectors y that appear in the 


interpolating helds. It is well known that a wide variety of interpolating fields 
can be used to represent a particular bound state within a reduction formula. We 
shall show in this subsection that the Bargmann-Wigner wave functions are the 
natural choice for these projectors in the heavy quark limit . 


As a hrst example consider a heavy vector meson. An interpolating held 0^1 
for a vector particle, with momentum P, should satisfy the condition 


X) 


(O|0/.(x)|P) = efxC 


—iP.x 


(90) 


where the polarisation vector satishes the transversality condition v.e = 0 and 
= - 1 . 

Thus, we can take as interpolating helds, for example, either 

0i(^) = NlT^lJg{x)-f^^pQ{x ), (91) 


1F) — 

= -N2j^T'iPq{x)Y'yi'ipQix) . 

Q 

The normalisation constants 


{0|ri/(,(0)|r-V-Q(0)|F) ’ 

and ^ 

" (opwoW*i(o)lA 


(92) 


(93) 

(94) 
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follow from the normalisation of the polarisation vector. 

Now one can use the equations of motion 

{ip — m)ip = 0 (95) 


for the quark helds to obtain 
if) — 

-J7^Tijg{x)Ylii^Q{x) 

Mq 

= [Ti^g{x)'y^ijQ{x) - 2iT-^q{x)—^ - ^q{x)] ■ (96) 

We see that the two interpolating helds differ by a term which can be neglected 
in the heavy quark limit. Thus in the heavy quark limit, with mg+niq = Mq, the 
two interpolating helds, 0^ and 0^, are equal along with the important relation 

Ni = N2 = N . (97) 

Note that the term neglected, Dj^/Mq, is exactly the kind of term which was 
dropped in deriving the HQET. One can also look at this result in another way. 

i'ipq{x)j^'ipQ{x) is in fact another possible candidate for the interpolating held 
for the vector meson but its overlap with the physical state becomes negligible 

_ ^_L 

in the heavy mass limit. Note also that i'ipq{x)^'ipq{x) corresponds to a p-wave 
contribution to the vector state coming from the anti-quark part of the heavy 
quark. In the zeroth order HQET recall that the antiquark part of the heavy 
quark held is suppressed and thus the overlap of this interpolating held with the 
physical heavy meson state becomes very small in this limit. 

The pseudoscalar case is simpler. One can show using equations of motion 
that the two interpolating helds 


0i(a;) = PiT0g(a;)750Q(a;), 

and 

if) — 

02(a;) = -P^^T'ipq{x)Yl5'ipq{x) , 
are equal along with the relation between the normalisation constants 

^ _ mq + TUq 

P2 Mq 


(98) 

(99) 


( 100 ) 


Note again that, as in the vector case, in the heavy quark limit, where mq+mq 

Mq, P1/P2 —^ 1 - 

Having established relationships between interpolating helds we now address 
the question: what are the consequences of the freedom of choice of interpolating 
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fields for a matrix elements involving heavy mesons and baryons? As an example 
we shall consider an arbitrary matrix element involving a heavy vector meson in 
the initial state 

M = {P’\J\P), (101) 

where the state \P) represents an incoming heavy vector meson, with momentum 
P and mass Mq and the hnal state is arbitrary as is the current J. Now use the 
reduction theorem to write this matrix element as 

M = hm - Mj) [ d^xe-^^-^P'\J(j)\x)\0 ), (102) 


The 0f(x) in (|102D where 0^ is given in eq. (j^) or (|^ . Thus one can write 


- Ml)N j, (103) 

J 


where xi^x) is either ^or and correspondigly N is either or 

The statement that one can use any one of the alternative interpolating fields 
means that they should give the same result in the reduction theorems. To see 

what this implies consider the second interpolating held, i.e x{dx) = the 

matrix element (|103|) . Integrating by parts we can shift the x derivative in y onto 
exp{—iP.x) to obtain 


M = lim (P2 _ mI,)N 2 [, (104) 

whereas with the other choice of the interpolating held we have instead 

M = d^xe-^^-^{P'\J{/)a^'ipQ{x)'ipqfi{x)\0). (105) 


Thus comparing we see that if the two expressions (|104|) and (|105|) are to be 
equal, then = A^i/A "2 as an operator inside the matrix element, at least in the 
mass shell limit for the heavy meson. In general this is not much of a restriction 
as the Ni are functions of the velocity v. However, we have shown above that, 
in the heavy quark limit, the ratio of the normalisation constants is unity and 
independant of the velocity, leading to ^ = 1. In fact we see from (|104|) and 
(|105|) that this condition means that only certain components of the quark helds 
entering in the interpolating helds contribute to the matrix element, i.e those 
satisfying 


Mx)]/ = iI)q{x) 

'j/'ipgix) = -ipq{x). (106) 


in the meson mass shell limit. In other words, in the rest frame of the meson 
on the mass shell, only the quark part of the heavy quark held and surprisingly 











also only the antiquark part of the light quark held contributes in the reduction 
formula. The same relation holds for the pseudoscalar particle. It is obvious that 
this physical requirement will be enforced by taking the interpolating held for the 
heavy vector meson, in the heavy quark limit, to be 

1 _ 7 /^ _ 

= (107) 

Q 

In fact even if we start with an arbitrary linear combination of the two helds, only 
the sum survives in the heavy quark limit. With this choice the matrix element 

picks up the correct projection operator to enforce the above condition ( |106|) . 

The above considerations turn out to be a general feature of the heavy quark 
limit. Given a particular interpolating held NT'ijjqTijjQ, then one can show, using 
the equations of motion, that this is equal to the interpolating held N^^T^jJq'y'^T^jJQ, 

upto terms of order 0{D^/Mq). Here, F is some Dirac matrix, possibly with 
derivatives. Thus the natural choice for the projection operator x{d) appearing 
in the interpolating helds for heavy mesons, eqs. (|69D and 0 ) , is 0 


Here, for example, F = 75 or F = /for the pseudoscalar and vector meson respec¬ 
tively. As we have noted earlier, the projection operator always appears in an 
integral over the space-time coordinates along with the relevant exponential con¬ 
taining the momentum of the physical state, corresponding to the interpolating 
held. Thus, by partial integration the diherentiation in the projection operator 
can be the shifted onto the exponential, so that we can always replace the x{d) 


of eq. (|I08D with 


X/3 


n = 


2jM, 


-[a + mp' 


(109) 


Q 


in all reduction formulae. 

Although we shown that the above form of the projection operators is the 
natural choice for the interpolating held for a heavy meson, we of course have the 
freedom to choose this form also for a light meson. This provides a uniform ap¬ 
proach for both light and heavy mesons . Similar considerations apply to the 


heavy and light baryon interpolating helds, containing the projection operators 

Xo-P')- 

®We have here introduced an explicit 1/yGTg factor in the projection operator to factor out 
the heavy mass scale. This is related to the fact that our states are normalised relativistically: 
{P\P') = 2E{2TTf5^{P - P'). 
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One now recognises that these projection operators (|109|) are nothing but the 
well-known Bargmann-Wigner wave functions. Hence, we have shown that the 
task of projecting out preferred particle states using products of spin half helds 
can be achieved in an elegant manner through the use of the so called Bargmann- 
Wigner wave functions. 

The Bargmann-Wigner wave functions are multi-spinor wave functions of 
some given rank and symmetry type, 


Xai...a„ {v) , 

which satisfy the free Dirac equation on all the labels 

(^- = 0 , 

.= 0 , 

(^- = 0 . 


( 110 ) 


( 111 ) 


The reason these wave functions describe particles of a given spin is amazingly 
simple [^]. In the rest frame the Dirac equation can be seen to be a projection 
by the operator (1 -|- jq). This means that the spin labels take on not four 
independent values but two. Put another way, this reduces SO (3, 1) down to one 
of its SU{2) factors. After imposing the Dirac equation on all the labels, the 
multi-spinor becomes a product representation of SU{2). On hxing to a given 
symmetry type one is actually considering an irreducible representation of SU{2) 
of dimension 2 s -|- 1 . 

For example, take x(n) to be totally symmetric of rank 2s satisfying il3). 
Then this is a wave function for a particle of spin s and even parity]^. In the rest 
frame, 70 is the parity operator for the spin 1/2 Dirac spinor, tensor products of 
it being the parity operator for the multi-spinor. 

A rank 2 symmetric bispinor, for example, describes a spin one particle as 
follows. The gamma matrices may be split into a symmetric class, and 
cr^jyC, and an antisymmetric class, C, 75(7 and The symmetric bispinor 

is then a linear combination 




( 112 ) 


The Bargman-Wigner equations ( 111 ) impose the constraints that 


= 0 , - v^<p^ , 


(113) 


so that, on taking cj)^ to be the polarisation tensor the wave function takes the 
form 

X{a/3} = ^ l^)^C]af3 = [P+^C]al3 ■ (H^) 

® Actually as defined only for 5'-waves. To get a handle on the P and higher waves allow the 
x{v, k) to depend on the orbital momentum 
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If one is interested, as we are, in using these multi-spinors to describe baryons 
and pseudoscalar mesons, then the above example is not enough. First, we may 
augment this with a rank 2 anti-symmetric bispinor X[ai 3 ]{v) describing a spin 
zero particle. Expanding this out in the three possible antisymmetric matrices 
and imposing ( p.ll|) we hnd that 

XM = ^[(1+ :^)75C']a/3- (115) 

Secondly, as the mesons are made out of quark anti-quark bound states the wave 
function has one covariant and one contravariant index. The contravariant index 
may be lowered with the help of the charge conjugation matrix C, so that one 
obtains a covariant bispinor whose symmetry or anti-symmetry may be declared 
unambiguously as above. We can then remove the charge conjugation matrix by 
multiplying by C~^ to obtain the projection operators for the pseudoscalar and 
vector mesons 


xJ- + 

(116) 

2^ 



(117) 




respectively, as expected. 

For the baryons our interest rests on the rank three trispinor Xa( 3 'y{v). After 
the Bargman-Wigner conditions are imposed each label effectively ranges over 
two values, so that this wave function has 2^ = 8 degrees of freedom. We wish 
to decompose this down to its irreducible parts (1/2 © 1/2 © 3/2), under SU{2). 
First set 

Xafiy — X[q:/3]7 T X{q:/3}7 ■ (H^) 

If we had not enforced the Bargman-Wigner conditions, then X[ai 3 ]'y would be 
reducible and the totally antisymmetric part would have to be projected out, 
that is, we would have to also impose 

X[o/9]7 T X[/37]a T X[7a]/3 0 • (H^) 

However, it is easy to see that after imposing the Bargman-Wigner equations 
this tracelessness condition is identically satished. X[a/ 3]7 represents a spin half 
particle. On the vother hand X{o/ 3}7 is reducible. To project away the totally 
symmetric part we impose the condition 

X{«/9}7 T X{l3'y}a T X{7Q!}/3 = 0 . (120) 

Let us keep denoting the rank three spinor of mixed symmetry that satishes this 
condition by X{ap}'i- It represents a spin 1/2 particle. We start with 6 degrees of 
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freedom from which we exclude the 4 coming from the totally symmetric compo¬ 
nent, leaving 2 components. The complete decomposition is then 


X.afi'y X[q!/ 3]7 T X{a/3}7 T X{a/37} • 


( 121 ) 


The reason for explaining this at length is that it allows for a rather nice 
description of the situation in the case of heavy baryons. Firstly, the totally 
symmetric component describes a spin 3/2 particle and so is suitable for the S* 
and n* baryons. Now for the heavy A and H baryons where the two light quarks 
are expected to form an antisymmetric diquark s = 0 state. If in (|121|) we let 
7 be the index on the heavy quark leg then X[ai 3 ]-y is just what we are looking 
for, and it has the form (following the same analysis as for the mesons), upto 
normalisation, 

A[q!/ 3]7 ■ ( 122 ) 

For the heavy S and hi baryons the light diquarks are expected to be in a 
symmetric s = 1 state. Still with 7 being the heavy quark label, the other 
spin 1/2 wave function X{a/ 3}7 i^i ( | 121 |) has the two light quarks in an s = 1 
combination. This then is the correct wave function for these baryons. There are 
two ways to get to an explicit form for this wave function. In the hrst method 
we follow the same route as we did for the mesons, namely we expand Xia/Sh 
terms of the symmetric gamma matrices in the a and (3 indices and impose the 
Bargmann-Wigner equations to hnd, again upto overall normalisation, 

A{a/3}7 = ^‘^7[(1 + , (123) 

with 

= O = (124) 

On imposing the tracelessness condition (|120|) , and after a little algebra, eq. 
(2.81) becomes 


A{a/3}7 


^[(1 + + ¥)ll.C]ap 

1[(7M + y/^)^5^x]^[(l + , 


(125) 


which is the wave function as put forward in |^, |^]. Alternatively we notice that 
a trace condition of the above form is easily solved in terms of antisymmetric 
objects, 

A{«/3}7 0[q: 7]0/3 T 0[^7]0a i (126) 


which on also demanding (| 111 |) is 


A{a/3}7 = 2 ^ 1^75C']a7W/3 + 2^^^ + f)lbC]0^Ua, ■ 


(127) 
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This is the form advanced in p3. The eqnivalence of the various forms of the 


wave functions was given in in a slightly roundabout way. The derivation 
here shows that from the Bargmann-Wigner (and group theoretic) point of view 
one form is as natural as the other. 

Likewise for the spin 3/2 wave function there are two natural forms. By 
expanding once more in the symmetric gamma matrices in the a and (3 labels, 
imposing the Bargmann-Wigner equations and then demanding total symmetry 
one hnds, up to normalisation, 


X { a / 37 } = <[(1 + ( 128 ) 

and also that ( p.24|) is satished with replacing 01). 

The second (and obvious) approach is to begin with the expansion in terms 
of the symmetric gamma matrices in any two of the labels and then add sym¬ 
metric permutations of these terms in all the labels. In this way we obtain, upto 
normalisation, 

X{a/37} = [(1 + 1^)7 mC']{«/?}< + [(1 + -f [(1 + ■ (129) 


3.4 Normalization to the lowest order 

Here we expand the interpolating field 0 q and the physical state $ in terms of 
powers of the inverse heavy quark mass, (pq and $ appear in the dehning relation 
for the normalization constant Nq, so that we are also led to an expansion of 
Nq. Some equations are spelled out in higher orders of (l/mq) to exhibit the 
expansion, and as we have need to refer to them at a later time, when we deal with 
the normalization at O(l/mo)- Using Eqs. (^) we can expand the interpolating 
helds 0Q, Eqs. (|6|) - (0), as 

0Q = 00 + X-01+ , (130) 

2mq 

where 0o is the interpolating field in the lowest approximation. To the lowest part 
00 we may associate a lowest order term in the mass expansion of the physical 
state. That is, we set 

I*) . + (131) 

From (1^ we have for pseudoscalar mesons (the extension to other mesons and 
baryons is immediate) 

A^Q(O|0g(O)x(n)0Q(O)|<h) = 1, (132) 

so that we are led to 
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where 


A'o(0fe(0)x(K)Q(0)|4»> = l, 


(134) 


and 


00(x) = 


Nl'il^q{,x)x{v)Q{x) 


(135) 


and similarly for the higher order terms. It is important to realise that to zeroth 
order Nq is independant of the flavonr of the heavy qnark. 

For a pseudoscalar meson we have shown above that we can take x{x) = 
^-^(1 + :j/) 75 . Take now v = (1,0) to get 


^(O|0,(O)7o75Q(O)|<h°) = 1 , (136) 

Mq 


where we have used the fact that upto zeroth order the heavy quark held is 
essentially Q+, i.e. it satishes 7oQ+ = Q+ in the rest frame of the heavy meson. 

We may relate the normalization back to the pseudoscalar decay constant fp, 
for we have 

{0\AomP) = fpMg , (137) 

so that 

(O|0,(O)7o75Q(O)|P°) = fpMg , (138) 


where fp is the zeroth order decay constant. Comparing with ( [136|) we have the 
result 

f°p = W ■ (139) 


N»qJMq 


Remembering that the zeroth order normalisation constant is independant of the 
havour, we have thus rederived the Voloshin- Shifman ^ scaling law for the 
pseudoscalar decay constant. 


4 Transitions from the Bethe-Salpeter Point of 
view 

In order to make use of dynamical models of bound state wave functions it is 
of practical advantage to cast the theory into a form which yields the above 
transition as an overlap of wave functions for the two heavy particle states. To 
do this, we hrst recall the dehnition of the Bethe-Salpeter wave function, which 
is the covariant generalization of the Schrodinger wave function. 

For simplicity let us concentrate on two-quark bound states for the derivation 
of the connection between current transitions and wave functions. The three- 
quark states are treated analogously. Consider a basis \P,a) for heavy mesons 
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consisting of eigenstates of the fonr-momentnm operator P, with a denoting the 
eigenvalnes of other commnting observables. The dehnition of the heavy-light 
Bethe-Salpeter wave fnnction ^^{P,xi,X 2 ) that we take is 

^^{P,Xi,X2)o,^ = {0\^|JQa{Xl)G{Xl,X2)^g{x2)\P,a) 

^^{P,Xi,X2)f^°' = {P,a\'iljQ{xi)G{xi,X2)iJqp{x2)\0) (140) 

where time-ordering is to be nnderstood here and snbseqnently. a,f3 are spinor 
labels and Q,q denote the heavy and light flavonrs. G{xi,X 2 ) and G{xi,X 2 ) 
are colonr matrices chosen to make the B-S amplitndes gange invariant. The 
canonical choice for them is the path-ordered exponential 

r^2 

Pexp[ig / A.dx]. (141) 

J X\ 

However we will not make an explicit choice yet. In onr case we are interested in 
the |P, a) corresponding to a 6 -meson or to a c-meson |<I)c) 

The Bethe-Salpeter wave fnnction is related to the two-body propagator |p4|P1 

/CQ(a;i,a; 2 ;a; 3 ,a; 4 ) = {Q\'ipQ{xi)G{xi,X2)'4)q{x2)'ipQ{x:i)G{x‘s„Xi)'4)q{xi)\Q) (142) 

by inserting a complete set of states between ipq and t/'q. In this way we get 

}Cq{Xi, X2] X3, X4) = <I)'5(P, Xi, X2)^^{P, X3, 0 : 4 ), {ti, t2 > ts, ^ 4 ), (143) 

IP,a) 


where ^^[P) and 4>'5(P) are the Bethe-Salpeter amplitndes dehned above. 

To get at the b ^ c transition, begin with the Green’s fnnction (c.f. (^)) 

G{xi,X2-,z]yi,y2) 

= {0\i>c{xi)G{xi,X2)i>q{x2)i>c{z).T.i>b{z)Myi)G{yi,y2)^/jq{y2)\0) 

= Y {0\'ipc{xi)G{xi,X2)'ipq{x2)\p,a){p,a\'ipciz).T.'ilJb{z)\q,b) ■ 

\p,a),\q,b) 

{(i,b\Myi)G{yi,y2)'Ggiy2)\0), (144) 


with xio > X 20 > Zq > yio > y 2 o and T any matrix (for {V — A) interactions 
it is 7^(1 — 75 )). Both snmmations are over complete sets of states. We wish 
to project ont the two-particle bonnd states. To do this, one overlaps with the 


all particle states that will appear will have only one given velocity, there is no chance 
of confusion in not including the corresponding momentum explicitly in the kets. 

^^From now on we drop the Dirac indices. We shall reinstate them whenever necessary. For 
example, Kq{xi,X 2 '■ X 3 ,X 2 ) in (1^) is a four spinor valued function There is no sum over 


Dirac indices in this equation. 
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appropriate Bethe-Salpeter wave functions and makes use of an orthonormality 
condition to obtain 


J d‘^xid‘^X2d'^yid‘^y2^p^{xi,X2)T{xi,X2;z;yi,y2)^P^{yi,y2), (145) 

where T{xi,X2', Z]yi,y2) is dehned by 
Gixi,X2;z;yi,y2) = 

J d^xsd'^x^d'^ysd'^yi ICcixi,X2; X3, x^)T (xs, X4; z- 1/3, |/ 4 )/Cf,(|/ 3 ,1/4; yi, 1/2), 

(146) 

so that T(x 3, 2:4; z; 1/3,1/4) is two-particle irreducible on both the (2:3, X4) and 
(1/3,1/4) legs. T(x 3, 2:4; z; 1/3,1/4) is sometimes determined in lowest order pertur¬ 
bation theory to be 


T{x 3 ,X 4 -,m,yi, •) = - x)S‘*\y 4 - z)S ■(x4,2/4) ® T, 


(147) 


where Sq{x 4 ,y 4 ) = (O|' 0 g(|/ 4 )V’g(a) 4 )|O) is the usual light quark propagator. For 
our purposes we keep T( 2 : 3 , 2 : 4 ; 2 ;; j/ 3 , j/ 4 ) in the implicit form (|146| ) for the present. 

As mentioned previously, the effective theory tells us nothing about conhne- 
ment, so that in particular it is inappropriate to use it in determining the Bethe- 
Salpeter wave functions (|140|) . It is in the region of interaction that the effective 
theory may be used. That is, it is reasonable to expect that the transition kernel 
T{xi,X 2 ] z]yi,y 2 ) may be accurately determined within a l/mg expansion for 
large mg and in a region where there are no hard gluons. 

T{xi,X 2 ', Z]yi,y 2 ) is to be calculated by using the dehning equations (|142 


1 144] , |146|) for the Greens functions and evaluating them within the effective held 
theory. Begin with the dehnition (144), replace the xjjb and tjjc by the zeroth order 


helds b and c and use (|47|) to go the uncoupled helds b and c to obtain 


g{xi,X2;z;yi,y2) 

= {0\^/Jc{xl)G{xl,X2)^/Jq{x2)i>c{z).T.^/Jb{z)^ljbiyl)G{yl,y2)i>qiy2)\0) 


= ( 0 |iy c{x,)G{xi,X2)ijq{x2yc{z).T.W 


1 -1 


w 


1 -1 


b{z)b{y4)W [III G{yuy2)My2)\0) 

= -{0\c{xl)c{z)\0){0\Tr'^jJg{x2)T^jJg{y2)G{xl,X2)W 
{0\b{z)b{yy\0), 


Xi z z yi 

Vc Vc Vf, Vf, 




(148) 


where in the last step we have used the fact that the transformed heavy quark 
helds, b and c, decouple from the gluons. The trace is over the colour indices. 
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Likewise we find for the two-body propagator, (|142| ) , 

lCb{Xi, X2] X3, X4) = (0l^b(Xl)G(Xi,X2)'lpg(x2)'lpb(x3)G(x3, Xq{x4)\Q) 

= -{Q\h{xi)h{x3)\Q){Q\Tr^q{x2)i!q{x4)G{xi,X2) 


w 


w 


n -1 


G(x 3 , 2 : 4 ) 10 ) , 


(149) 


where there is, again, a trace over the colonr indices. By inserting these into the 
expression (|146|) , one is able to solve for T(xi, X 2 ; |/i, 2 / 2 ) 

T(x3,X4;z;i/3,y4) = - J d'^Sid'^S2S'^{x3 - z) 5 '^{y 3 - z)S~^{vc; X4, Si; Xi, X3) 

g{xi, Si; 2 / 1 , S2)S~^{vb, S 2 ,2/4; 2/3,2/i) , (150) 


where 


SqiVc] X2,X4;Xi,X3) 

= {0\Tr'lpq{x2)'ipq{x4)G{Xi, X2)W 


XI 


w 


X3 


-1 


G'(x3,X4)|0) , (151) 


and 


<Sg(vb;y4,y2;d3,yi) 

= (0jTr^q(y4)^g(y2)G(y3, y4)W 


ys 

J 


W 


yi 

L^i> J 


n -1 


G(|/i,|/2)|0), (152) 


g(xi,X2-,z-,yi,y2) 

= (0jTr'i/^g(x2)i’g(y2)G(xi, X2)W 


XI z z y\ 

Vc Vc Vb Vfj 


^( 2 / 1 , 2 / 2 ) 10 ) (g)r 


(153) 


Again in the above equations the traces are over the colour indices. The S ^ are 
dehned through 


I d^X4Sqivc;x2,X4;x,,X3)Sq\v,;x4,s;x4,X3) = 5\s-X2) (154) 

and 

J d^yiSq\vb; s, 2 / 4 ; y 3 , yi)Sg{vb] 2 / 4 , 2 / 2 ; 2/3,2/i) = <5^(s - 2 / 2 ) • (155) 

An important result follows from this. We see that in the heavy quark effective 
theory one can factorize the Green’s function and so solve for T, a situation that 
is not available in the full theory, except at the perturbative level. 
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4.1 The Form of the Bethe-Salpeter Wave Function 


A comparison of (150) with (^) or (|86|) allows us to infer the Bethe-Salpeter 
wave function of the heavy meson to be of the form 

$®^(a;i,X2) = Xo2\v)A^sixi,X2), 

(156) 

or in momentum space 


^'if‘(P,k) = xJ(v)Al{P,k), 

(157) 


for some multi-spinor function A. This may be obtained directly from the LSZ 
representation of as follows. Writing the state \P,a) in (|14(]|) in terms of the 
interpolating field (^), one sees that the B-S amplitude for a heavy meson can 
be written as 


^a^{Xl,X2) 

= Xp\v)N lini (p2 - M‘^){Q\ijQa{xi)G{xi,X2)^'^{x2) 

j dV*'^'^'0Q(2/)V’95(l/)|O) 

= Xp\v)A'f^{xi,X2). (158) 


The form ( [156| ) follows when we take the heavy quark fields, to be the 
corresponding fields in the effective theoryf^. Thus when we transform to the 
decoupled effective fields (c.f. 13) we get 


4^(11,12) 


. j dSe-^P-yQ\y)W 



^ 9 . 5 ( 2 /) 10). 


X\ 

V 


G{xi,X2)i:^ix2) 


(159) 


Finally the ^ in the decoupled heavy quark propagator is set to unity (i.e. 6^) 
by the ^{1 + ]/) projector of the y, using the same procedure as in deriving (|85D , 
leading to the form of the B-S amplitude given in (|156|) and (|157|) . 

The B-S amplitude (|156D (or |157D has the property that the heavy quark 
label a satisfies the Bargmann-Wigner equation. (See footnote 4). This single 
important fact leads to enormous simplification of matrix elements. 

follows from the discussion in the subsection on the Barg mann-Wigner wave functions 
and interpolating fields that the form of the B-S amplitude (1^) is valid also for light mesons, 


if we take Q to be a light quark. It only simplifies on going to the heavy decoupled effective 
fields. 


38 
















4.2 Universal Form Factors and Symmetries of the Effec¬ 
tive Theory 


Returning once more to the transition (pq), we can write it in the following 
compact form (it makes no difference at this point if we use (|85D rather than (^) 
for the general features we are about to derive) 


V ,), (160) 

with 

E{vc,Vb) = 4:NbNc j (fkd^qM{q, k). (161) 

Here r(r) is either 75 or ^ (75 or f') depending on whether we are dealing with 
a pseudoscalar or vector particle in the initial (hnal) state. Also, because of the 
heavy flavour symmetry present in the zeroth order HQET, the normalisation 
constants are equal upto 0{1/M) viz. Nb = = N. We also note that H is 

independent of the heavy meson masses as we have pulled out the mass scales 
and therefore it is independent of flavour. 

The most general form for H is 


•^{vb,Vc) = Fo{w) + + l^bF2{w) + , 


(162) 


with oj = Vb-Vc- However, because of the projection operators and , the 
matrices that go to make up 5 may as well be replaced by the identity so that 
(|160|) takes the universal form 


MbM,i{w)Trf{ 




(163) 


with i{w) = Fq — Fi — F 2 + F 3 . This is a very strong result, for it says that 
regardless of whether the transition is pseudoscalar to pseudoscalar, pseudoscalar 
to vector or vector to vector, there is only one form factor. The differences 
between these various possibilities for transitions rests completely in the different 
r’s and r’s that appear in the trace. This is what is referred to as the spin 
symmetry present in the heavy quark effective theory. Thus doing the traces we 
have 


{D{vc)\cjxb\B{vb)) 

{D*{vc)\c-fxb\B{vb)) 

{D*{vc)\c-fx'y5b\B{vb)) 

{D*{vc)\c'yxb\B*{vb)) 

{D*{vc)\c'yx'y5b\B*{vb)) 


MBMD^{w){Vb + Uc)a 

^JMBMD*^{w)[{l + w)el - Vb.e*Vcx\ 

^Mb*Mb* i{w){{vb + Vc)xe*.e - Vc.eel - Vb.e*ex] 

^ Mb* Mb* i{w)ie*''{vb + VcYe^^e^p^x (164) 
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Though the physics of this situation is clear, i.e. spin information is at least 
of order 0(l/mQ), let us also predict this spin symmetry result from a different 
more formal approach. Consider once more the transition as expressed by (|6^ , 
and vary the c quark field according to the transformations (^51). As the action 
is invariant under this change, one obtains a Ward identity 

($'= ($,|^,75/Morv^„|<f)b), (165) 


where the state is the one that corresponds asymptotically to ys/MoX- On 
taking to be the polarisation vector for a spin one particle with four velocity 
u, Mq = 1 and y to represent a pseudoscalar, we see that the new state is a 
vector complete with the standard form of the wave function! This allows us to 
relate the different y factors, but the reader may be wondering about the overall 
normalization Applying the spin symmetry argument to ( |136|) , shows that 
the normalizations Np, Ny are the same at 0(1). 

In this case then, ( |165|) relates the pseudoscalar to vector transition with the 
pseudoscalar to pseudoscalar transition, and the universal nature of heavy quark 
transitions is manifest. Also the Bethe-Salpeter wave functions themselves can be 
shown to have a universal multispinor factor A. For example, the vector meson 
Bethe-Salpeter wave function <l>y is related to that of the pseudoscalar meson <I)p 
via spin symmetry: 


(0|75^Q(a:i)hF 


Xi 

V 


G{Xl,X2)^|Jq{x2)\P) 

75 


{ 0 \Q{xi)w 


Xi 

V 


$ 


Q 
V 5 


G(Xi,X2)^g(x2)iV} , 
(166) 


where \P) and \V) denote pseudoscalar and vector meson states, respectively. 


4.3 Normalization at 


2 

max 


One immediate consequence of the heavy quark effective theory is the normal¬ 
ization of the single form factors, ^{w), in the matrix elements (|164|) , at or 
w = 1. This normalization is achieved by comparing with the conserved flavour 
charge number[^ 


(*q(^)I / <i‘xi'Q(x)-yoi’Q(x)\<iQ(P)} = 1, 


(167) 


where Q denotes a single heavy flavour. This equation is exact in the sense 
that it is valid beyond the effective theory. Thus the zero component of the 
vector part of the transition Q —> Q at zero recoil is absolutely normalized to 
any order in (l/mg). This normalization is also valid if we take the limit that 
rriQ —>• oo, where we can replace the heavy fields ^Jq^x) by their effective fields 

this equation and in the rest of this subsection, the 1 on the right hand side stands 
symbolically for 2E{2Tr)^S^{0) 
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Q{x). Now we prove, using the flavour symmetry of the action Eq.(^), that this 
normalization at q^ax extends also to heavy flavour changing transitions at 0 ( 1 ). 
Let the unitary matrix U of Eq. (^) be 


U = 






and the heavy quark in Eq. (|167|) be the charmed quark c. Thus we consider the 
transformation c —*• + 5c, c —*• + 5c. From LSZ reduction we learn that we 
can replace also the c and c in the states by the unitarily transformed ones to get 

1 = j d?x{ph{x)+5c{x))'^o{l3h{x) + 5c{x))\^jib+Sc) 

= J d^xb{x)'yob{x)\^b) + {55y{^c\ J d^xc{x)'yoc{x)\^c) 

+ {f3(5){55) Jd^xb{x)'yoc{x)\^c) + (*hc| Jd^xc{x)'yob{x)\^b)^ ■ 

(168) 


As U is unitary, + (55) = 1. Using this and the fact that = {Mb — McY 
is below threshold, so that the amplitude (<hc| / d5’xc{x)'')Qb{x)\^b) is strictly real, 
we conclude that also the zero component of the vector part of the transition 
5 —>• c is absolutely normalized at q^ax 0 ( 1 ), that is 


($,(P)| j d^xc{x)^oh{x)\^b{P)) = d 


(169) 


Comparing with the first of the equations (|164|) leads to the wellknown result, 
e(l) = l, toO(l) 

As the action (|3lD is invariant under the transformations (^5[) , we have also a 
manifest spin symmetry of heavy quark transitions, as exhibited in ( 1 ^) at 0 ( 1 ). 
Using this we obtain absolute normalizations of spin related transitions. For 
example, the pseudoscalar to vector transition is normalized by taking T = ^570 
with e = (0, 0, 0,1) and v = (1,0) in (|165|) . Thus, we have 


{V,{P)\ I d^xc{xhs^,b{x)\Pb{P)) 

= {P,{P)\ J d^xc{xhob{x)\Pb{P)) = 1 (170) 


using the fact that 70 acts on (1 + 7o)/2 to give (1 + 7o)/2. 

In Section 5 we will derive the remarkable fact that this normalization holds 
also to 0 (l/mQ). 


5 Effective theory to higher orders 

In this Section we will have a look at some of the interesting features and con¬ 
sequences of the effective theory at 0(l/mQ). As the derivation of the higher 
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order form of the effective theory is no more difficult than the derivation of the 
hrst order form we will see firstly how to derive all the higher order terms, in the 
effective theory. The 0{l/mQ) term is determined explicitly. Our aim then, in 
the rest of this Section, is to apply the formalism we have just derived to hrst 
order in the heavy quark mass. 


5.1 Derivation of Higher Order Action 

Let us write the action in terms of the helds Q (or Q) as 

OO 

So = S» («) + x: S5(«)/(2mo)” 


(171) 


n=l 


where Sq{v) was given in (^) and derived from the standard action in (^9]). In 
order to determine the Sq{v) we need to substitute (^7]) into (|56D and keep terms 
up to 0(1/(2771^)"^). To simplify the formulae we dehne 


P^ = p — 'j/v.D 


= YDt, 

(172) 

which has the useful property that 


1- 

o 

(173) 


is the covariant derivative transverse to the velocity. 

The action (^) is then 

= + (174) 

On setting Sq{vq) = QOnQ, from (|174|) we can hnd the On to be 


a,= 




k\ 




(175) 


Exercise: Check this equation. 


This is rather too implicit. As an example we now show how to put Oi in a 
conventional form, 

= + (176) 
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where 


(177) 


= |l7p.7j and [D'‘,D‘'] = -igF‘“'. 

The action up to hrst order is then 

Sq = + (178) 


Returning to the question of why does the Foldy-Wouthuysen transformation 
have the form that it does, we note a tell-tale feature already present at 0(l/mQ). 
This is that there are no derivatives in Oi in the v direction. This is important for 
if there were such a derivative, there then would arise terms in diagrams, where 
this operator was inserted, that would go like v.P = tuq in the numerator thus 
raising this 0{l/mQ) correction to 0(1). Quite generally the appearance of 
in the transformations ensures that no derivatives in the direction of v appear at 
any order in the 1/mQ expansion. One can see this quite easily from the equation 
(|174|) . There the only term containing such a derivative is 


^i{P^)/{2mQ)-^^jj^i{P^)/{2mQ) 
We now use the identity 


(179) 


e»yle-« = A+-Y. JB, |B, -[B, Tl]]. 


(180) 


n=l 


n times 


with B = i{p^)/2mQ and A = iv.D. In this case [B,A] = . Thus 

to all orders one does not get derivatives in the v.D direction. 

The operators On are appropriate for the Q helds. Equivalent operators On 
may be found for the Q helds from 


S^{vq) = QOnQ = QOnQ. 


The required relationship is 


On{x) 


w 


'On(x)W 



(181) 


(182) 


and as the On only involve covariant derivatives, this relationship may also be 
expressed as 

dn(A) = On(A') , (183) 

where the gauge held A' is given in (^9]) and satishes v.A' = 0. 

In order to get an action which completely disentangles the high and low 
frequency states (i.e. Q+ from QJ) requires another set of transformations. We 
will not need that action and leave this as an exercise to the interested reader 
(see 0). 
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5.2 Normalization of the interpolating field at 0{l/mQ) 

From the normalization of the interpolating field (|132|) , the expansion of Nq (|133|) 
and of the physical state (|131|) , and taking the action correction Si into account, 
we get for mesons 


= 


Q 




Q, 


2m, 


Q 


Si 




Q 


Here Nq can be read off to be 


Nh = -iK? 




(184) 


-Q + iJ.xQSil^^) + (Ol^xQI®!)) . (185) 


If there are H-corrections (see Sec. 3.2), they must also be incorporated in 


5.3 — corrections to form factors 

mq 

In this subsection we consider the 0(l/mQ) corrections to the zeroth order meson 
transition form factor ^{w). Here we will establish that there are no 0(l/mQ) 
corrections to the normalization, at zero recoil, of the form factor that we found in 
Section 4.2. The strategy is the following. One relates, once more, the vector part 
of the transition at equal velocity to the c-number charge. As this has already 
been normalized to 1 at 0(1), any corrections at 0(l/mc) necessarily vanish. This, 
together with the spin symmetry, yields the desired result. We give two proofs 
of this result. The first is purely diagrammatic, while the second is rather more 
explicit. In the second method we explicitly calculate the 0(l/mc) corrections 
to the mesonic matrix elements of the 5 —> c vector and axial transition currents 
away from the zero recoil point and show that they vanish at zero recoil. The 
0(1/mb) correction can be evaluated in an analogous way. 

There are three possible sources of corrections that need to be taken into 
account. One stems from the extra term that appears at this order in the action, 
namely Sq. The second correction is due to the fact that the helds i/'c that 
appear in the reduction formulae are not just the lowest order helds c but rather 
(1 +iP'jr/{2mc))c. The last contribution can come from perturbative corrections 
to the H-functions. 

From the methods of proof it will be apparent that both l/ruc and 1/mb cor¬ 
rections vanish at threshold. We therefore only address the 0(l/mc) corrections 
directly. 


5.3.1 Diagrammatic Proof 


This proof was given in and elaborated upon in Let all of the corrections 
be thought of as insertions into the lowest order diagram and represent these 
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insertions by a cross. The cross stands for either a momentnm kick or a sonrce 
of glne or both. For the 6 —> c transition this is exhibited in £g. 2. Note that the 
correction at the vertex is only on the c-qnark side. Also there shonld be (dressed) 
gluon lines running between all the quark hnes[^ which for ease of visualization 
are not drawn. The corresponding charm charge normalization diagrams that 
should sum to zero are exhibited in £g. 3. There are two crosses on the vertex 
corresponding to the corrections on either side of the vertex. 

The charm charge normalization diagrams may be expressed as twice those 
of £g. 2 where once more at the vertex the cross is a correction on the right 
hand side only. This doubling up occurs because we may read the diagrams 
in any way we choose as the momentum flow is the same on either side and 
rr( 7 ^j .. .. . 7 ^ 1 ). But these last diagrams are just those of the 

zero component of the vector part of the b ^ c transition. Hence that component 
of the transition must vanish. To show that the axial component also vanishes, 
one invokes the fact that the spin symmetry also implies a normalization condition 
(see Section 2.7). 

For this proof to work, we need to make the assumption that there is a sensible 
I/ttiq expansion of the H-function (see Sec. 3.2). We have no control over this 
term, and it would require some ‘modelling’ to ascertain its form. This lies outside 
the scope of these lectures and we take as a working assumption, that such an 
expansion exists. 


5.3.2 — correction to vector and axial transition matrix elements 


211. Let us hrst look at the 


In this subsection we follow the proof given by Luke 
correction coming from the fact that, upto 0{l/mc, in the reduction formulae 
we do not have just the held c but (1 + /{2mc))c. This hrstly leads to a 

modihcation in the currents 


cTh —>• cTh — 


2m, 


■ci pj_ Tb, 


(186) 


where the T subscript is with respect to Vc- Secondly the held corrections will 
also enter in the interpolating helds (p.30| ). These correspond to corrections to the 
physical state ( |131|) . The total held corrections to the yccF (when one transforms 
to the free helds) portion of the LSZ reduction formula are thus 


^ (^x{vc)P^{x)c{x)c{z)T - x{vc)c{x)c{z) p^ {z)T^ , (187) 

the hrst term being an interpolating held correction, the second is the correction 
to the current (vertex correction) already written in eq. (|18(j| ). The charm prop¬ 
agator is made up of p and the identity matrix, while xi^c) projects onto the 

the gauge vt-A — 0 there are no gluons at all connecting the heavy quarks. In a general 
gauge these do contribute, but always in the form of Wilson lines. 
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(1 + component. One may then move the to the right hand side of the 
hrst term without engendering any new gamma matrix structure. So as far as the 
gamma matrix structure is concerned both terms are of the same type. One can 
do better, and establish that the wave function correction vanishes completely 
when one goes to the mass pole. This is not needed for our subsequent analy¬ 
sis; but the reader might like to keep in mind that for all intents and purposes, 
1$^) = 0 in (|131|) . Thus we will only consider the vertex correction, i.e. the 
correction due to the modihed currents. Let us now calculate this for the B ^ D 


and B - 
element 




transitions. Thus we begin by considering the following matrix 
1 


2m, 




(188) 


where T^ = 7^(1 — 75). To evaluate this consider hrst the matrix element 


{D{v,)\{aDx r^6)(0)|S(ufe)). 


(189) 


From our previous analysis, in Sections 3 and 4, we see that, because D\ does not 
contain any gamma matrices, this matrix element can be written immediately as 


^MsM^Tr75(^^)r^(^^)75 [Av,^ + Bv,^ + 

= ^jMBMoivb + Vc)^,{Avbx -h Bvcx) 

-C[gf,x{l -w) + VcbiVbx + VcxVbbi] - CAie^j^xpaV^v^ , (190) 


where A, 5, C are unknown fuctions of tc = Vb-Vc- Now multiply by and use 
the equation of motion 

nVc- D= —mcC]^^ (191) 

to write 

- mc{D{vc)\c]/j:'^b\B{vb)) = ^JMbMd{Aw + B - C){vb + Vc)p . (192) 

On the left hand side we now notice that will be hit by the projector on 
going to the free (c) held, and therefore can be set to unity. Recall the discussion 
at the end of Section 3.1 where the projector has already been used to remove 
the from the c-quark propagator. Hence we get 


- mc{D{vc)\cT^b\B{vb)) = \JMbMd{Aw + B - C){vb + Vc)p . (193) 

Comparing with the hrst of equations (|164|), we get 


Aw + B — C = —md . 


(194) 
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We also have, upto Oil/mq), 

{D{vc)\idx{B: yh)\B{vb)) 

= {Pb-Pc)x\cT^h\B{vb)) 

= [Mev^x-MDV c\){vbPVc)^i{w). (195) 

The left hand side of this equation can now be written as 

{D{vc)\{a Dx T^b)\B{vb)) + {D{vc)\{cr^i D b)\B{vb)) • (196) 

The hrst term in the above expression we have already evaluated in eq. (|19CI|) . 
Hence multiplying eq. (|195|) by and using the equations of motion ivb-Db = 
rrib'^b^ 3.nd the same argument we used to derive eq. (p.93|) we arrive at 

A + Bw — C = {Mb — mb)i{w) — Mdw^{w) . (197) 

Comparing with eq. (|194|), we get 


(198) 


A{1 + w) — C = A^{w), 
A-B = 

where A = {Mb — rub) = {Mr, — rric). 

Now let us hnally look at the vertex correction 


1 


2m, 


■{D{vc)\ci Tp\B{vb)) 

{D{vc)\ci Dx 'yiTp\B{vb)) 

^MbM.Tn, ^ ^ ^,{Avbx + Bv,x + C^x) • 


2mc 

1 

2mr 


(199) 


Here the functions A, B and C are the same as in eq. (|190|) . We see immediately 
that only the vector transition current contributes and the trace can be worked 
out to give 

- -^\/MbM„{vb-Vc)f,{A{l + w)-3C) , (200) 


which, on using eq. 


, becomes 


2m, 


■,jMhMc(vb - Vc),, (2.4(1 + lo) - 3A^(u')) 


( 201 ) 


The second kind of correction comes from the fact that we have to include 
the effects of the new terms which appear in the action at order 0(l)/mc 

1 


5i = 




( 202 ) 
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Inserting this in the path integral leads to a correction equal to 


^^{D{v,)\ld-'xT(c(iD^fc) (i)(cn6)(0)|B(«i)) 
+±-SD(v,)\ISxT (i)(cr„t)(0)|B(t.j)) (203) 


Using the reduction formulae and going to the free (tilde) helds, as in Sec. 3 we 
hnd, since {iD^Y does not include 7 matrices, that the hrst term will lead to the 
trace 

= y/MbMc^^{vb + Vc)ixVi.w) , (204) 


where ri{w) is an unknown function. This simple form appears because the 
appearing in the two free c-quark propagators can be put to unity because of the 
projector x(nc) appearing in the reduction formula. 

In the second term we follow the same procedure but we see immediately 
that because of the the projector x{vc) arising from the reduction cannot act 
directly on the second of the free c-quark propagators. However we recall (p^) 
that near the heavy quark mass pole the main contribution in the propagator 
comes from 

„• + i^c) 

Vc-Pc - rric ' 

Thus the second term leads to the following trace 





* ~ (1 + (1 + (1 + A 


-0'1/A- 


75^^ 


z/A 


(206) 


It is easy to see that the only independant form of is 


, 


(207) 


leading to the correction 

[D{w -1) + 3E]{vb + Vc)fx, (208) 

where D and E are unknown functions of w. 

Thus the total 0{l/m^ correction to the B ^ D decay is 



(7 -|- D{w 


1) -h 3E) {vb + Vc)^ + (2H(1 + w)- 3A^) {vb - n^)^ 


(209) 

We can follow the same procedure for the other s-wave matrix elements to 
obtain hnally the following set of matrix elements for flavor changing currents 
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upto Oil/rric). 


(L>(vc)|c7/,6|5(v6)) 


\fMbMc (^^{w) + - l)-D(w) + 3E{w)}^ {vb + Vc)a 


+ ^;^{2(1 + w)A{w) - 3A^{w)}{vb - Vc)^ 

{D*{vc)\c-ff,b\B{vb)) 

A 

2m^ 

{D*{vc)\c'-i^')^h\B{vb)) 


= ^MbM^ 


1 

2mr 


(1 + + ■T^i'niw) - E{w)) ie*X<e!^P<TM > 


(210) 


( 211 ) 


= ^MbM, 


<UH (1 + w^) 


2mc J 2mc 


-Vb.e'Vc^ |?(w)(l + “ 2/l(ui) + D(w) - E(w)) 


+t>i.e*JV^^(2-4(m) + D(wy) 


( 212 ) 

(213) 


Here e*‘' is the polarisation vector of the E* meson. 

To determine the normalisation at zero recoil it is necessary to consider the 
neutral current matrix elements, which are governed by the same set of form 
factors because of the heavy flavour symmetry. Following the same procedure as 
for the charged currents, we get 


2 1 
^(w) + + (^^^ - l)T)(t(;) + 3E(w)) (vi + V2)^ 


(D(v2)jcj^cjE)(vi)} 

= ' 2m, 

(D*(v2)ic^f,cjE>(vi)} 

= M, [(1 + A^)«u>) 

^11 Lq 


(214) 


( 2 ? 7 (t(;) + 2(1 + t(;)A(t(;) - (1 - t(;)Zl(t(;) - 2 i?(t(;)) ie*''viV2e^p^f^ 


2m^ 

(L>*(n2)|c7^75c|T>(ni)) 


(215) 


= M, 
1 


^(w) ((1 + w) -4(1 -m;)^ 

V 2mc/ 

(^2(1 — w‘^)A{w) — (1 — w’^)D{w) + 2(1 + w)E{w) + 2(1 + w)r]{w)'^ 

( / /(A \ 1 


-Vi.e*V2p. I^(tf) + 


4A 

2mr 


+ 


2m, 


■ {2r]{w) — 2(2 + w)A{w) + wD{w) + 2E{w)) 
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+ vi.(^*V2i,^^(2A(w) + D(w)) 
{£>*(t> 2 )|^„c|C’(t'i)> 


Me 

[(^(w^) + 2 ^ (^(^) (0 + V2)f,ei.e2 


1 

2mc 

{2A{w) - D{w)) (ui -f V 2 )biVi.e 2 V 2 .ei 


^(tc) - (2?7(tc) + 2(1 -I- w)A{w) -I- (1 — w)Diw) — 2Eiw)) 

1 2mc J 

^*2^,V2■el 

^(tc) - {2r]{w) + 2(1 -I- w)A{w) -l- (1 — w)D{w) — 2E{w)) 

1 2mc 

> eif,vi.e*2 


( 217 ) 


Here vi and V 2 are the velocities of the incoming and ontgoing D-mesons, respec¬ 
tively and the ei and 62 are the corresponding polarisation vectors. 

We observe immediately that all the 0(l/mc) corrections vanish in both the 
flavonr changing and nentral matrix elements at the zero recoil point w = 1, 
except those terms proportional to rf(w) and E{w). Next we recall that the 
nentral vector cnrrent is normalised at tc = 1 since at this point, ni = ^2 = n, it 
is a symmetry current corresponding to the conserved c-quark number. Thus 

{D{v)\c'yoc\D{v)) = 2McVo = 2McVo + 3-^(1))^ (218) 

and 

{D*{v)\c'yoc\D%v)) = 2M,vo = 2M,vo (l + ^(2h(l) - 2^(1))) . (219) 


Here we have used the fact that ^(1) = 1. Thus one gets the normalisation 
condition 


rj{l) = E{l) = 0. 


( 220 ) 


Thus we arrive at the conclusion that all corrections to the zeroth order result 
vanish at the symmetry point w = 1. 

One can calculate the baryon transition matrix elements upto 0(l/mc), in 
a similar manner, using the baryon Bethe-Salpeter amplitudes developed ear¬ 
lier. Without going into the details we give the results, upto 0(l/mc), for the 
transition A;, —> Ac |@, 1 ^, JT^I, HU, HU, HU. 


{Ac{vc)\c-f^b\Ab{vb)) 

= u{vc)-fi,u{vb) I^a(m^) + ^ (hA(w) + A^a(u;))| 



1 -1- w 


u{Vc)u{Vb) 


( 221 ) 
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and 


{Ac{vc)\cj^'y5b\Ab{vb)) 

= HvchfiAMvb) lu{w) + (va{w) w^)A^a(w) 


1 A^a(w) 
rUc 1 + w 


2mr 


Vbf,u{Vc)j5u{Vb) . 


1 + w 


( 222 ) 


Here, as in the mesonic case, ^a{w) and r]A{w) are unknown functions normalised 
as follows, at the zero recoil point: 


a(i) = 1, 

hA(l) = 0. 


(223) 


6 Renormalisation and the Relationship with 
QCD 

So far we have considered how QCD is related to the HQET theory at the classical 
level. That is we have hxed the tree level coefficients in ( |171|) by relating the 
theories by a Foldy-Wouthuysen transformation. When it comes to renormalising 
the theories a discrepency arises. One obvious difference is that the higher order 
operators that appear on the right hand side of (p.71|) are non-renormalisable. Due 


to inhnities that the insertion of such operators entails the tree level coefficients 
will all need to be renormalised. Indeed as the operators in the expansion have 
increasing mass dimension, the higher the order in l/mq one goes to the worse the 
divergence becomes. In order to renormalise therefore, one will need to perform 
an inhnite number of ‘experiments’. That problem aside, for the moment, even 
the HQETq differs from QCD once loops are taken into account. 

It might come as a surprise that this is so, as we have derived the HQET from 
QCD, and one might have expected that knowing the renormalisation constants 
of QCD would be enough. The differences arise because the transformations used 
do not respect the regularisation that is needed to dehne QCD in the hrst place. 
To see how this situation comes about consider the integral 




(p2 


m 


2 '|2 
Q) 


(224) 


which, when regularised with dimensional regularisation, apart from a pole has 
a logarithmic behaviour hnrtq/p. Now let us perform the integral by expanding 
about the mass shell p = mqv + k, 


<rk 




4:mq{v.kY \ rnqv.k 


+ ... = 0 . 


(225) 
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(We have used the fact that / d'^kk'^ = 0 in dimensional regularisation). The 
discrepancy lies in the fact that, while at tree level we can keep the k (and gluonic 
momenta) bounded by mg, we are unable to do so in loops. Similar manipulations 
show that even for convergent integrals one gets a mismatch. Notice that the type 
of terms that appear in the series expansion are of the form of the higher order 
operators that appear in the HQET. 

Is this the death knell for the HQET? No not at all. This is a common situation 
in effective held theories. The naive effective held theory will correctly reproduce 
the physics at some scale (in our case the long-distance physics, k -C mg) but 
will require corrections to reproduce the physics at other scales (here the short 
distance physics or for typical momenta k > mq). Indeed one can give a system¬ 
atic treatment of the short-distance corrections. As the ehective theory is QCD 
in disguise, the ‘experiments’ that are needed to hx the renormalisation constants 
in front of every term on the right hand side of (|171j) are simply a comparison of 
the ehective theory and QCD up to the given order at some scale where we will 
demand equality between the two theories. We will come to terms with what this 
means presently, but hrst we wish to establish, for external momenta k -C mq, 
that at leading order QCD and HQETq are indeed simply related. 


6.1 Renormalisability of the HQET 


0 


The lowest order ehective theory, HQETq, has a power counting renormalisable 
action (essentially as the counting is the same as for QCD). In principle we can 
stop here and conclude the theory is renormalisable. There is a quicker way to 
arrive at this, namely that, in the gauge v.A = 0, the heavy quark decouples 
altogether from the rest of the theory. The heavy quark is free, while the rest of 
the theory is renormalisable in this gauge. 

If the reader does not particularly like this gauge, then he or she can work 
in a covariant gauge, again with free heavy quarks but with the inclusion of the 
Wilson lines (see section |2.5| ). Wilson lines are basic as observables in QCD and 
better be renormalisable. Indeed this was postulated for such operators of smooth 
loops by Polyakov and proved shortly after by and . A review of these 
matters may be found in |4^ . 


6.1.1 Wavefunction Renormalisation of the heavy quark 

Let us calculate the wavefunction renormalisation in the Feynman gauge, as we 
will have need of it later on. From now on quantities with asterix superscripts 
are defined in HQETq. The same objects without the asterix are the equivalents 
in QCD. We wish to evaluate (with p = mv -|- k) 


(1 + ¥) 
2 




_ I _ 

+ q) — mq^ 
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(1 + ^) 4 ,2 . f 1 

2 3^ ^ J {2n)^ v.{k + q)q'^ 


One now makes use of the identity, 


—- = / da 
a^b 


a 


n—1 


to rewrite the self energy as 

JllAify. 

2 3^ ^ , 


da 


0 {aa + 

d^q 1 


(27r)"^ {av.{q + k) + q^Y 


(1 + ^) 4 


4 % 

3 f/*V ( 4 ^) 2 -,/ 2 r(^/^) /q da{a^/A-av.k) 


-e/2 


In order to perform the last integral one scales 

a —Av.ka 

and then changes variables to 


= 


3 a 


to arrive at 


(1 + ^) 4 




2 3 " (4 

= (^^(-4..A:)^-T(1 - e/ 2 )r(e - 1 ). 

The singular part is easily extracted from this expression. It is 

2 3 167r2 e 

from which we see that the wave function renormalisation is 

_ 16 //*Vl 

The renormalised inverse propagator is, therefore. 


^(2,0) “ * 


.(1 + 7 ) 


v.k 1 


16 g 


*2 


3 167r2 


In /i + ... . 


(226) 

(227) 

(228) 

(229) 

(230) 

(231) 

(232) 

(233) 

(234) 


Now the anomalous dimension for the heavy quark 7 q can be determined from 
the renormalisation group equation 




^(2,0) “ 0 • 


(235) 
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The beta function is of order 
obtain 


and so we can ignore it for present purposes to 


'^Q 


8 9 *^ 

3 167r2 ■ 


(236) 


Before closing this section let us return to the observation that what we are 
really calculating are expectation values of Wilson loops. The heavy quark wave- 
function renormalisation should be the composite operator renormalisation of the 
Wilson loops and consequently the anomalous dimension ( P36|) ought to be the 
anomalous dimension of the Wilson loops with end points. Indeed it is as has 
been shown by |^. The reader should bear this older literature in mind. Many 
of the calculations performed in the context of HQETq could be extracted from 
past work on Wilson loop renormalisation. 


Exercise: Prove (|227|) and hll in the details of the calculation leading to 


Exercise: What of the renormalisation of the 




component of the propagator? 


6.2 Comparison of the HQETq and QCD 

In the previous section we calculated the one loop anomalous dimension of the 
heavy quark in the lowest order HQETq. On the other hand, the one loop anoma¬ 
lous dimension for a quark in QCD is, 


4 9 ^ 

3 167r2 ■ 


(237) 


So there is a mismatch between the QCD and the HQETq predictions. While 
renormalisability for the HQETq is not an issue, it is its relationship with QCD 
which requires some elaboration. 

In the sequel, for simplicity, we will assume that the theory of interest is that 
of the heavy quark coupled to the glue, with no other matter present. The general 
situation can be dealt with along similar lines but is rather more involved. 


6.2.1 Comparison of Gluonic IPI Diagrams 

One of the important features of the HQET at lowest order is that there are no 
heavy quark loops at all. This is easy to see. Heavy quark loops correspond to 
the fermionic determinant, which is gauge invariant, and therefore we can work in 
the gauge v.A = 0 to conclude that the determinant is gauge held independent. 
Consequently there are no heavy fermion loops. An alternative way of saying this 
is that, in the rest frame, the heavy quark can only propagate forward in time, 
leading again to no loops. So how does this square with QCD? 

In fact what is relevant here is a very important result due to Appelquist and 
Carazzone |Q(see also [^). Consider an n point function with only external 
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gluon legs whose momenta are of the order of A; {k « mo). These authors tell 
us that, in renormalised QCD, for such diagrams, those which contain a heavy 
quark loop are supressed by at least a factor of k/mq relative to those which do 
not contain such a loop. The detailed proof of this theorem for a simpler model 
plus a useful guide to the literature may be found in the book . 


Convergent Diagrams 


In order to get a feel for this part of the Appelquist Carazonne theorem, 
consider a diagram with n > 5 where all the external gluon lines are connected 
to the heavy fermion loop and where all sub-diagrams are convergent. One has 
the generic expression (we drop 7 -matrices, traces and so on). 


i=l i=l i=l qi 

F 1 


Adi 


)(n 


i=i P - E;=i qj - rriQ p + ELi kj - tuq 


-), (238) 


where r^(gj) is the gluonic IPI F-point function, (Fig.4). 

Since the degree of divergence of T'^ is 4 — F, and as the subdiagrams are 
all to be convergent, we require F > 5 as well. Scale the integration variables 
in this formula as qi —>■ mqqi and p —>■ mqp and note that, as everything in 
sight is convergent, the scaling law of T is T^imqqi) = [q^). The scaling 

argument shows us that the diagram behaves as 


CoTUq "■ CikrriQ 


(239) 


for some hnite constants cj. Such an n-point function without a heavy quark 
loop would, by dimensional analysis, go like 


ck^-^ 


(240) 


validating the claim. The essence of the scaling argument that we have used is 
simply that, because of the presence of the mass in the fermionic loop, we can 
set the external momenta to zero without generating a divergence. 

At the other extreme, consider the situation where none of the external lines 
land on the heavy quark loop, (Fig.5). This time we wish to get a handle on 


n Ai)<5(^F)F”+^(fc,gi)(II 

i=i qi 


Adi 


2=1 


2=1 


i=i P - ELi qj - mq 


)• (241) 
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Suppose that 5 < F < n. In this case the heavy quark loop behaves like 
as long as the loop momenta q are small. In any case contract the loop to a 
point and consider the reduced diagram. The reduced diagram’s overall degree of 
divergence is {F — n) < 0, so that it superhcially converges. All the subdiagrams 
converged and so we may conclude (by Weinberg’s theorem) that the reduced 
graph converges. The shrinking of the loop to a point is the mq oo limit 
of the original diagram, so that we hnd that such graphs behave as 
which are once more supressed with respect to diagrams without fermion loops. 
When {F — n) > 0, one estimates this integral by looking at the regions where 
all the integration momenta are small (relative to mq) which is the dominant 
contribution (again by application of Weinbergs theorem). One can show that 
such diagrams behave as These are clearly supressed. 


Divergent Diagrams 


We now have to address the question of power counting divergent diagrams. 
According to Appelquist and Carazzone, once one renormalises the basic diver¬ 
gent graphs then their inclusion into other diagrams is subdominant. The argu¬ 
ment for this is that, heuristically, the subtracted diagrams (at the scale /i mq) 
behave as (here A is a cutoff) 


/i 


[In — + 0{—)] - [In A + o( JL)] = o{ — 

mq mq mq mq mq mq 


(242) 


An example of this behaviour (in dimensional regularisation) is afforded by the 
one loop correction to the gluon self energy due to a heavy quark loop. This is 


= {k^ku - k‘^ri^t,)T{{k) 


(243) 


where 


U{k) 


2 9^ 
3167r2 


- 7 — 6 / dxx{l 

€ Jo 


x) In 


mq + k'^xil 
Anfj? 



(244) 


Once we subtract Yiik"^) — n(/i^) and substitute into diagrams, the contribu¬ 
tions are bound to be convergent (this is the essence of the renormalisation pro¬ 
gramme), as the high k behaviour is tempered. On the other hand, for k mq, 
the (subtracted) logarithm will behave like 


9 ’ 


(245) 
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as in ( |242| ) . 

Once these, and similarly subtracted gluonic three and four point functions, 
are considered in the above expressions (with F < 4), one obtains the estimates 


X 0(— ) or 0(— ), 
mg mq 


(246) 


again establishing the suppression of such heavy quark contributions. 


The Relationship between QCD and HQETq 


We have thus shown that, in a particular (mass dependent) subtraction scheme, 
the n-point light particle IPI diagrams r„( 5 f, mg, m,/i) in QCD and the n-point 
light particle IPI diagrams P* ((?*, m*, /i) in the HQETq are related by the follow¬ 
ing simple relationship which holds to order l/mq and for /i -C mq, 

Tnig, mq, m, /i, k) = m*, p, k). (247) 

(Here we have reintroduced all of the other light, relative to the heavy quark in 
question, particles). 

The couplings which appear on the right hand side of this equation are those 
that would be obtained in QCD with the heavy quark omitted. One immediate 
and very important consequence of all of this is that the beta function of QCD 
at scales below the heavy quark mass goes over to the beta function of QCD 
without the heavy quark. This is the idea behind grand unihed theories. 

There is a technical aside that should be made and that is that Appelquist 
and Carazzone work in the Landau gauge. The reason for this is that in any 
other covariant gauge the gauge parameter itself needs to be renormalised. The 
theorem has been carried over to other covariant gauges. We have already seen 
that the gauge v.A = 0 can, in this respect, be problematic. 


Cautionary Remark: It is worth emphasising that one should not blindly apply 
the decoupling theorem. Here is a counterexample. Consider the divergence of 
the axial current at one loop. The anomaly is mass independent and 

does not vanish no matter how heavy the heavy quark is. 

When choosing the renormalistion prescription we do not wish to introduce 
any extra mq dependence into the renormalised couplings. It is thus good policy 
to employ some MS scheme for which the renormalisation constants are mass in¬ 
dependent. We can relate the mq dependent scheme of Appelquist and Carazzone 
to an MS scheme in the following way. Let the unrenormalised, but dimensionaly 
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regularised, n-point light particle IPI diagrams be rf/^„(5fo, rngo,e)- There 
relationship to the renormalised F, mgj, m,,/i), determined in two schemes 
i = 1,2, (one might like to think of i = 1 as the mass dependent subtraction 
scheme used by Appelquist and Carazzone, while i = 2 is some MS scheme) is 

Tu,n{9o,mQo,mo,fi,e) = zf^Ti^n{gi,mQi,mi, fi). (248) 

By construction, the rj^„( 5 fj, mg*, m*,/x) are hnite and are related by the hnite 
ratio 

T2,ni92,mQ2,m2,lJ.)=\^—j ri,„(5(i,mQi,mi,/i). (249) 

Similar formulae hold for the HQETq n-point functions as well 

/ \ 

(^j (250) 

Now starting from (|247|) one can pass to an MS' scheme on both sides to obtain 
(we have switched notation so that the MS scheme couplings are to be understood 
and we hope our laxness causes no difficulties), 

^n{9, rriQ, m, p, k) = Z{g, tuq, , m*, p, k) , (251) 

for a certain function Z{g,mQ, fi) which is given in terms of the hnite renormal¬ 
isation constants Zi and Z*. 


Exercise: What is the relationship between Z{g,mQ, g) and the renormalisation 
constants Zi and Z*7 


6.2.2 Comparison of Green Functions with 2 Heavy quark legs 


One would hope that a simple relationship like (|251|) holds also for external 
heavy quark legs and indeed it is true. Here, however, we should specify that 
the momentum running through the heavy quark line is “almost” on shell. That 
is we should take that momentum to be mqv^ + k^, where k -C mg. This will 
ensure that we can pass to the HQETq (we need this condition even at tree 
level). Given the discusion above, we will take it for granted that there is no 
need to take into account heavy quark loops and that the renormalisation of the 
light part of the theory has been accomplished. Again we should distinguish the 
case of convergent versus divergent diagrams. Also, given that the HQETq is 
renormalisable, it is apparent that we only need to renormalise the heavy quark 
propagator and the heavy quark-heavy quark-gluon vertex, as we must do also 
in QCD. The following argument, which gives us our sought for relationship is a 
variant of that due to Feinberg and to Grinstein []^ . 
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Convergent Diagrams 


Consider within QCD, Green functions G 2 ,ni with two external heavy quark 
lines. Denote by G 2 „ the same Green functions evaluated in HQETq. To simplify 
life later on we take it for granted that there is a projection of (1 + y)/2 on these 
both in QGD and in HQETq. Set n > 2. The overall degree of divergence of 
such Green functions is (1 — n) < 0, so that, apart from subdiagrams which are 
divergent, these Green functions would be convergent. Gonsider any convergent 
diagram contributing to such a G{ 2 ,n)- 

According to |^, the integrals of such diagrams are dominated by the region 
where the loop momenta satisfy qi < K, where K <C mg is the typical momentum 
flowing through any of the external gluon lines or k. If this is the case, then one 
can substitute all the heavy quark QGD propagators with the corresponding 
HQETq ones, 

1 1 . ^ 

^-^ -T-, -^-• 252 

f+q-rriQ + q) - niQ 

The two diagrams will then agree up to sub dominant contributions from wich we 
may conclude that for such convergent diagrams G( 2 ,n) = ^{ 2 ,n)- 

As an example of the argument presented in p| consider the diagram £g. 6. 
We wish to determine the behaviour of 


i^ — _ I _ 


(253) 


Split the region of integration into a region where is less than and one 
greater than so that symbolically we have 





(254) 


In the integral /q^ the momentum q is small relative to mg and may be ignored 
in the QGD quark propagator. Furthermore, up to subdominant terms, the QGD 
propagator becomes the HQETq propagator. So this part of the integration 
coincides in the two theories. 

It is easy to see that the integral d'^q, on the otherhand, is itself subdom¬ 


inant for s < 1 and smg > K. In this region one can set A; = fci = 0 in ( |253|) and 
change variables to g —>• mgq, to find 


noo 1 11 roo 

/ ^ 

IsmQ fmg + q- mg g^ mg Js 


W A 


'+4-1' q 


(255) 


This is clearly subdominant and shows us that the integration momenta which 
he somewhat above K are irrelevant. 
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Divergent Diagrams 


What of G 2 ,n for n = 0,1? Focus on the one loop correction to the three 
point functions r( 2 ,i) and r* 2 i)- Within dimensional regularisation the integrals 
are hnite and so up to Oil/mo) they agree. However, both integrals diverge 
as we approach four dimensions and need counterterms. Unfortunately there is 
no guarantee that the counterterms preserve the relationship. If we differentiate 
with respect to the residual momentum or the external gluon momentum, the 
differentiated diagrams are hnite. Then the previous arguments imply that 

d d 

^r(2,i) = —1(2,1)*+0(/c/mQ,g/mQ) 

d d 

^r(2,i) = —Ti2,iY+0{k/mQ,q/mQ). (256) 

So, apart from counterterms, the two r( 2 ,i) are in agreement. Now the countert¬ 
erms are of the form ar °2 i) and a*r( 2 '’^^, where the superscript 0 indicates the 
tree diagram, and a and a* are inhnite constants. Now one chooses a and a* to 
ensure equality (up to subdominant terms). However, such a choice is bound to 
be niQ dependent. So at this point we have 

rg,.) = r;«, (257) 

where R denotes the renormalised Green functions. 


The Relationship Including Two Heavy Quark Lines 


We can pass from this mq dependent renormalization scheme to some mass 
independent scheme (as we saw above) by multiplying the Green functions by 
some hnite renormalisation constant, so that 


= C{mQ/Y,gs)Z{mQ,Y,gs)Tl^i){,k,q]Y) 


* R 


(258) 


After this long discourse we hnally arrive at the generalisation of (|251|) that we 
sought; 


T( 2 ,n){g,mQ,m,Y,k) = G(mQ//x, ^)Z(^, mg,/i)'^/^r( 2 _„)(^*, m*, p, fc). (259) 

The reader may be surprised that we have by-passed a discusion of the heavy 
quark two point function altogether. The reason for this is that, in deriving 
(|259D, we only need to know about the wave function renormalisation to get the 
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multiplicative factors. These are fixed between the two and three point functions 
by gauge invariance. Notice also that we have an equation for C{mQ/fi, g). That 
is, by applying the renormalisation group both to the left and the right of (|259|) , 
with 77 , = 0 , (and ignoring beta function contributions) we obtain 

-2(70 - 7yc 

(260) 


dC 


6.3 Matching 


The HQET is constructed to reproduce correctly the low energy behaviour of 
QCD. That is, in any diagram where the external momenta are much smaller 
than the mass of the heavy quark one can pass from QCD to HQET, remembering 
to multiply by the relevant hnite renormalisation constants. This means, in 
particular, that we can hope that the mismatch should not depend at all on very 
low energy physics, such as infrared singularities, physical cuts etc. However, 
as we pass to lower and lower energies perturbation theory becomes more and 
more unreliable and the equality of QCD and HQET more and more difficult to 
ascertain. For example there may be non-perturbative effects that ‘see’ the quark 
masses. These questions aside, it is apparent that there is a range of energies for 
which the HQET and QCD will agree to the desired accuracy in inverse powers 
of mg. 

Mismatches between QCD and HQET certainly begin to set in as one ap¬ 
proaches scales of mg, for in this regime the decoupling theorems are no longer 
applicable. In order to get agreement between QCD and the HQET one ‘matches’ 
them at /i = mg. These matching conditions manifested themselves in (|251|) , 
where the extra multiplicative renormalisation constants arose to take into ac¬ 
count hard gluon exchange, missed by the decoupling theorems. One extends this 
to other operators. This means that if we have some operator A and calculate in 
QCD or in HQETq then the relationship is 


(H(mg)) = C'o(mg,/i)(Ho(/7))o + 


Ci(mg,/i) 


2 m, 


Q 


(24i(/i))o -|- ... . (261) 


The notation is that on the righthand side one is working within QCD, while on 
the left hand side all expectation values are taken with respect to HQETq and 
the subscript j on the operators indicates that they are the operators together 
with the held and action corrections at order (mg)“T The Wilson co-efficients 
Cj{mQ,fi) are dehned by this relation. 

The idea that one should match operators in this way goes back to Voloshin 
and Shifman H].See also E^. The following example is the one originally worked 


out in 
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6.3.1 Matching for qVQ 

Let us see that the Wilson coefficients are calculable in practice. We write 


{qV'ipq) = Cl{mQ,^){qTQ)Q +... . (262) 


The matrix T can be either 7 ^ or 7 /^ 75 . We wish to get a handle on Cq g* 


Act on both sides of (p62|) with d/d/i, to obtain 


/i 


^^0 

-r- = wr 

djji 




0 ) 


(263) 


where 7 r and 7 ^ are the anomalous dimensions of the currents qT'ipq and qVQ 
respectively. 

In QCD the operators qT'ipq, with massless quarks and T = 7 ^, are conserved 
currents, and are partially conserved currents when T = 7 ^ 75 . They have van¬ 
ishing anomalous dimensions. The currents correspond to flavour symmetry in 
QCD. Adding mass terms for the quarks breaks the symmetry softly. This means 
that the ultraviolet properties of the theory are not changed. In loop integrals 
the presence of masses certainly effects the soft momentum region of integration, 
but they play no role in the ultraviolet region. Consequently, the currents remain 
finite and at most require a finite renormalisation. Thus, even if the quarks have 
mass, the anomalous dimension is zero. We are left with 


= (264) 

The reason that there may, nevertheless, be a non-zero anomalous dimension 
for the effective theory operator qVQ is that in the HQET it does not correspond 
to a (partially) conserved current. One can see this by noting that the heavy 
quark kinetic term is totally different to the usual light quark kinetic term and 
so there is no hope of a symmetry at this point. Consequently it is not protected 
from renormalisation. In practice such a renormalisation is required, as in QCD 
the currents exhibit a logarithmic dependance on mq which diverges as mq —> 00 , 
that is it diverges as we pass to the HQET. 

The solution to ( ^64|) has been discussed in other lectures at this school and 
is 

= exp (- [ ^ ^d^'^^^Co^(/io,^*(/io)), (265) 

\ P[9 ) ) 

where g is the running coupling defined by 


and with initial condition 


, dg{p') 

^ dg,' 


fd{g{p')) 


9[p) = 9* ■ 


(266) 

(267) 
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Take /iq = mQ in (|265|) to give 


CoimQ/^i,g*{fi)) = exp 



In order to determine the Wilson coefficient we need to know two things, 


7 r(5(') and 

Co • 


z 


(269) 


zz 


Both may be calculated perturbatively provided that /i and mg are large enough 
to ensure that ^(/i) and gimg) are small (i.e. A -C /r -C zjiq). 


{i) The Anomalous Dimension yf 


We can determine the anomalous dimension of the current by demanding that 
the renormalisation group equation 



(270) 


holds. We also know from (|236D and ( |237| ) that 


4 g^ 


(271) 



3 16ir2 


and therefore we only need to determine the pole part of Tp to work out 7 p. The 
diagram of interest, £g. 7, corresponds to 



(272) 


One can set /c = 0 without engendering an infrared divergence, and as we are 
only interested in the pole behaviour, we do so. The in the numerator gives a 
convergent integral and therefore we ignore it. Furthermore, the pole part of the 
integral must be proportional to and thus it is given by 



(273) 


We have then 



(274) 
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( 275 ) 


Plugging this into (|270|) (and once more ignoring the (3 term) we obtain 


-,*2 


7r = 4 


1671 ^ 


{a) The Wilson Coefficient Cg 


At tree level, QCD and the HQET agree. So we have 

Co ( 1 , 9{mQ)) = 1 + 0{g{mQf), 


(276) 


which will suffice for our needs. We only want to work 
logarithms. 

Recall that 

m = s (-'>•4^+ (lij)) 

with 6o = 11 — where n/ is the number of ‘active’ quarks. Below the 6-quark, 
but above the c-quark, threshold we have Uf = A etc.. Putting all the threads 
together we can determine the exponent in (^65|) 


at the level of leading 


(277) 


We thus arrive at 


where 


and 



4 rSM dg' 

bo Jgimq) g' 



9{mQ) ■ 


^0 


exp 



girng)) 


V «^(/^) ) 



2 


a = — 


bo' 


(278) 


(279) 


(280) 

(281) 
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If one wants more precision one can go to the next to leading logarithms 
approximation. For this one would need to calculate 7 p to the next order. 


Exercise: Determine C from (|260|) . 
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A Symmetries of the Heavy Quark Action 

In this appendix we give the most general symmetries of the heavy quark ac¬ 
tion.There are 16 symmetries, 8 “vector” and 8 “axial vector” types. In order 
to exhibit this we need to introduce some notation. Suppose the dimension of 
space-time is d. Pick a 7 matrix algebra for this space and let the 7 matrices 
be n X n matrices. As a basis for the 7 choose ^ and Any n x n may be 

expanded in terms of a basis P generated by the 7. The basis we choose is such 

that P = P+ (g) P- where 

{^,P+} = 0 (282) 

and 

[^,P_]= 0, (283) 

depending solely on whether there are an odd or an even number of y-*- in the 
product used to dehne the given element of P. Clearly P+ and P_ each have 
dimension n^/2. The symmetries of QySv.dQ are 

SQ = (a+P+ + a-p_)Q 

SQ = (5(a+P+- a-p_) (284) 

Notice that ^ the P+ symmetries are “chiral”, though, obviously, none of them 
is anomalous. 
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B Figure Captions 

Figure 1 Born diagram for the electron-proton potential 

Figure 2 Insertion of 0(l/mc) corrections in a 6 —> c transition 

Figure 3 Insertion of 0(l/mc) corrections in a c —> c transition 

Figure 4 Gluon n-point function where all the external gluon lines are connected 
to a fermion loop 

Figure 5 Gluon n-point function where none of the external gluon lines land on 
the heavy quark loop 

Figure 6 Gonvergent four vertex integral 

Figure 7 1-loop heavy light vertex correction 
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